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I.  ELECTRICAL  AND  OPTICAL  PROPERTIES 

HOT  CARRIER  DISTRIBUTION  FUNCTION: 

EXPERIMENTAL  DETERMINATION,  DESCRIPTION,  AND  APPLICATION. 
W.  Pinson  and  R.  Bray 


I .  Introduction 

The  amount  of  IR  radiation  transmitted  through  p-type  Ge  can  be  modulated  by 
applying  a  strong  electric  field  or  by  changing  the  temperature  of  the  Ge.  In  both 
cases,  the  changes  in  transmission  are  attributable  to  a  change  in  distribution 
function,  hence  a  redistribution  of  the  holes  in  the  upper  two  valence  bands.  In 
the  following  both  effects  are  combined  to  yield  the  energy  distribution  function  of 
the  hot  carriers.  The  distribution  functions  have  then  been  applied,  with  the  aid 
of  the  low  field  transport  parameters  of  the  Ge,  to  compute  a)  the  energy  loss  rates 
to  optical  and  acoustic  phonons,  b)  the  hot  carrier  mobility,  and  c)  a  parameter 
describing  the  impact  ionization  rate  for  deep  acceptors. 

I I .  Energy  distribution  function  of  Hot  Carriers 

Fig.  1  shows  the  relevant  parts  of  the  band  structure  of  Ge  and  the  expression 

for  the  interband  free  carrier  absorption  as  given  by  the  theory  of  Kahn  and  Kane^ . 

The  IR  absorption  Is  due  e.g.  to  a  photon  of  energy  hv  causing  a  hole  at  energy 

in  band  1,  to  jump  vertically  to  band  3,  Equation  (1)  on  the  figure  shows  that  the 
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absorption  coefficient,  a,  is  proportional  to  the  square  of  the  matrix  element,  H ^ 
for  transitions  from  bands  1  to  3,  to  the  hole  density,  p  at  e^  and  inversely 
proportional  to  hv  and  to  f,  a  function  of  the  curvature  of  the  valence  bands.  The 
absorption  coefficient  may  be  modulated  by  changing  p  («^),  either  by  applying  a 
field  to  the  Ge  or  by  changing  the  temperature.  If  the  other  elements  in  equation 
(1)  are  Independent  of  field  and  temperature,  then  the  absorption  coefficient  at  a 
given  X,  depends  only  on  the  distribution  function  of  the  holes.  For  the  measure¬ 
ment  of  absorption  as  a  function  of  temperature,  the  temperature  of  the  crystal 
changes  slowly,  the  holes  are  always  In  thermal  equilibrium  with  the  lattice  and  thus 
have  a  Maxwellian  distribution.  When  a  pulse  of  high  field  is  applied  to  the  material 
the  holes  attain  a  steady  state  distribution  but  are  not  In  thermal  equi 1 ibrium  wi th 
the  lattice.  As  we  shall  see  this  steady  state  distribution  is  no  longer  Maxwellian. 

Fig.  2  is  a  block  diagram  of  the  apparatus  used  to  measure  the  change  in 
transmission  with  field.  IR  from  the  monochromator  is  focussed  on  the  Ge  and  then 
refocussed  on  the  detector.  A  5  psecond  long,  high  field  pulse  applied  to  the  Ge 
causes  the  carriers  to  redistribute  themselves  In  the  upper  valence  bands  (1  and  2) 
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and  thus  modulates  the  IR  transmitted.  The  redistribution  time  is  much  faster  than 
the  rise  time  of  the  voltage  pulse.  The  modulation  shows  up  as  a  pulse  on  the 
detector  and  scope.  The  normally  transmitted  intensity  "ly"  is  measured  in  the 
same  system  by  chopping  the  IR  at  13  cps.  All  field  measurements  were  made  with  Ge 
held  at  77°K. 

Fig.  3  shows  the  percentage  modulation,  HA  ly/iy  for  fields  of  130  and  2800 
volts/cm  the  lowest  and  highest  used.  A  lT  is  the  change  in  1^.  caused  by  the  field. 

The  negative  values  are  due  to  an  increase  of  hole  density  far  out  in  the  heavy  hole 
band;  the  positive  values  due  to  a  decrease  close  to  k  *  0.  The  ranges  corresponding 
to  1  «  3  and  2  •*  3  transitions  are  marked  on  the  figure.  Intensity  modulation  approach 
ing  50^  was  observed.  We  have  found  that  the  percentage  modulation  also  depends  on 
the  crystal  orientation,  on  the  direction  in  which  the  IR  traverses  the  Ge  and  on  the 
direction  of  polarization  of  the  IR  relative  to  the  field  direction. 

Fig.  4  shows  the  change  in  absorption  coefficient  calculated  from  the  percentage 
modulation  for  a  field  of  800  volts/cm.  For  comparison,  we  also  show  in  this  figure 
the  change  in  absorption  at  zero  field  for  a  temperature  increase  of  100°K  above  the 
base  temperature  at  77°K.  Comparison  of  the  two  curves  shows  qualitative  but  not 
quantitative  agreement  between  the  absorption  coefficient  changes  due  to  temperature 
and  field.  This  we  attribute  to  the  non-Maxwel l ian  character  of  the  field-perturbed 
distribution.  Had  the  field  simply  heated  the  carriers  to  a  higher  "temperature** , 
then  the  field  curve  would  coincide  with  some  temperature  curve.  No  such  correspon¬ 
dence  could  be  found!  The  non -Max wel l ian  energy  distribution  of  hot  carriers  may  be 

obtained  by  properly  correlating  the  absorption  measurements  a(T.  )  made  with  varying 

2  L  ^ 
temperature  ,  and  those  made  by  applying  electric  field  pulses,  a(E),  Use  is  made  of 

the  theory  that  the  absorption  coefficient  at  a  given  wavelength  is  proportional  to 

the  hole  population  at  a  specific  energy  in  band  l ,  irrespective  of  whether  the  hole 

population  is  determined  by  varying  the  lattice  temperature,  or  by  the  applied  field. 

A  simple  proportion,  as  given  in  equation  2,  then  permits  determination  of  the  hole 

populations  when  the  field  is  applied.  Thus, 

px<f)  -  <» 

where  the  hole  population  at  any  specific  energy  and  known  crystal  tenperature, 
p^CT.)  is  given  by  the  Maxwellian  distribution  of  Eq,  3  in  the  absence  of  applied 
field. 


(3) 
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where  p  Is  carrier  concentration,  <x  is  absorption  coefficient,  the  subscript  X  Is 
the  wavelength  for  which  the  particular  quantities  are  determined,  and  E  and  refer 
to  electric  field  and  temperature  measurements  respectively,  it  remains  to  specify 
what  energy  In  band  1  corresponds  to  the  particular  wavelength  used.  From  the 
proportionality  between  p  and  a,  and  using  the  Maxwellian  relation  of  Eq.  3,  we 
obtain 


‘W  TL/2  ~  e"thh/kT|- 


(4) 


3/2  2 

A  semi  log  plot  of  a(T^)  vs  1/T^  taken  at  selected  wavelengths  yields  the  cor¬ 

responding  energy  relations  between  in  band  1,  and  X.  An  attempt  at  a  similar 
analysis  for  band  2  was  not  successful,  in  that  no  distinguishing  straight  line  slopes 
could  be  obtained  at  any  X, 

Fig,  5  shows  some  calculated  distribution  functions  for  field  strengths  between 
130  and  2150  volt/cm.  The  energy  range  between  the  dashed  lines  is  that  for  which 
there  is  reliable  data,  i ,e.  for  which  application  of  Eq,  4  yields  reasonably  straight 
line  slopes,  hence  values  of  The  energy  distribution  for  800  volt/cm  is  shown 

again  In  Fig.  6,  where  it  is  conpared  with  a  Maxwellian  distribution  at  169°K,  which 
has  the  same  average  energy  as  the  800  volt/cm  distribution.  The  areas  under  the  hot 
carrier,  and  the  Maxwellian  distribution  curves  are  the  same  within  experimental 

error.  It  Is  to  be  noted  that  the  average  energy  is  surprisingly  low  for  such  a  large 

3 

field  strength.  The  non-Maxwel 1  Ian  character  of  the  hot  carriers  is  demonstrated  , 
Coopered  to  a  Maxwellian  distribution,  the  hot  carriers  show  a  depletion  of  population 
at  low  and  very  high  energies,  and  a  piling  up  of  carriers  at  intermediate  energies, 
where  the  concentration  becomes  greater  than  can  be  found  in  a  Maxwellian  distribution 
at  any  tenperature.  The  depletion  at  high  energies  is  attributed  to  the  large  proba¬ 
bility  of  emission  of  optical  phonons  of  energy  k9  *  .037  ev  by  those  hot  carriers 

opt , 

with  *  >  kfi  . 

opt 

ill.  Application  9f  JMstrlbutlao  Function?. 

A.  Energy  Loss  Rate  Calculation. 

Two  qualitative  features  stand  out  from  the  above  results;  (1)  the  non-Maxwe! I  Ian 
character  of  the  distribution  of  hot  carriers,  and  (2)  the  low  average  energy  even  at 
strong  electric  fields.  The  latter  result  implies  the  presence  of  a  very  strong 
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energy  loss  mechanism  which  can  only  be  attributed  to  strong  interaction  of  holes 
with  optical  phonons.  The  deduction  was  quantitatively  tested  as  follows:  The  energy 
loss  rate  for  optical  and  acoustical  phonon  interactions  was  calculated  utilizing  the 
experimentally  determined  distribution  function,  represented  by  p(e)  in  ec\iat!on  5: 

(if).  “  J*  <  (*?)  -(*?)>  p(,)  d*  (5) 

loss  o  \  '  em  \  /  abs 

where  the  two  terms  in  brackets  represent  the  difference  between  energy  loss  and  gain 
rates  in  emission  and  absorption  processes;  Ac,  the  amount  of  energy  lost  (or  gained) 
by  the  hole  in  an  interaction,  is  k®opt  for  optical  phonon  interactions,  and  fita  for 
acoustical  phonon  interactions.  The  former  Is  a  constant,  while  the  latter  has  to  be 
taken  as  a  function  of  scattering  angle  and  carrier  energy.  The  brackets  indicate 
averaging  over  scattering  angle  where  necessary.  The  collision  time  t,  different  for 
the  emission  and  absorption  processes,  requires  knowledge  of  the  appropriate  deforma¬ 
tion  potentials  for  optical  and  acoustical  phonon  scattering.  These  were  taken  from 
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previous  analysis  of  the  low  field  mobility  .  The  calculated  energy  loss  rate  may  be 
conpared  with  the  experimentally  measured  one.  The  latter  is  equal  to  the  energy 
gain  rate  of  the  carrier  in  steady  state,  which  is  Just  given  by 

(df/dt)ga!n  “  «  * 

where  p  Is  the  measured  mobility  at  given  field,  f:,  and  p  is  the  total  hole  concen¬ 
tration. 

Conparlson  of  the  measured  and  calculated  results  at  various  field  strengths  are 
shown  In  Fig.  7.  The  solid  line  Is  the  experimental  result.  The  energy  loss  rates 
celculated  are  indicated  separately  for  optical  and  acoustical  phonon  interactions, 
it  is  seen  that  the  calculated  loss  rate  to  optical  phonons  is  very  close  to  the 
measured  total  loss  rate,  and  two  orders  of  magnitude  larger  than  the  loss  to  acousti¬ 
cal  phonons,  which  may  indeed  be  neglected.  The  agreement  is  best  where  the  distri¬ 
bution  function  is  determined  with  the  best  precision,  at  about  800  volts/cm.  The 
change  in  absorption  was  somewhat  different  for  the  sample  oriented  in  the  (100)  and 
the  (111)  directions,  so  separate  distribution  functions  were  calculated  for  the  two 
directions  and  the  theoretical  energy  loss  calculated  for  the  resulting  distributions. 
The  power  loss  Is  slightly  higher  for  the  (111)  than  the  (100)  direction. 
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The  optical  phonon  Interaction  had  previously  been  introduced  to  explain  the 
temperature  dependence  of  the  mobility  of  holes,  but  the  strength  of  the  interaction 
had  been  empirically  determined  to  fit  the  data.  The  same  parameter  now  enables  us  to 
account  quantitatively  for  the  carrier  power  loss  and  the  relatively  low  average 
energy  of  the  hot  hole  distributions. 

B.  Impact  Ionization  Rate  of  Deep  Lying  Acceptors 

The  distribution  curves  show  relatively  few  carriers  at  energies  much  above 
Thus  It  should  be  difficult  to  ionize  impurities  with  deep-lying  levels,  «|. 
This  has  already  been  observed  experimentally^,  We  can  look  at  this  more  quantita¬ 
tively  as  follows.  We  assume  that  the  ionization  cross  section  per  carrier,  a,  Is 
proportional  to  the  difference  between  the  ionization  threshold  and  the  carrier  energy. 
The  Ionization  rate  is  also  proportional  to  the  carrier  velocity,  v  J7.  If  the 
product  a  v  «<(•  -  f|)  Is  now  averaged  over  the  distribution  function  p(t)  from 
«l  to  •  (see  equation  on  Fig,  8)  the  result  Is  a  value  which  Is  proportional  to 
the  impact  ionization  rate.  The  same  product,  O  v,  may  be  averaged  over  Maxwellian 
distributions  at  several  temperatures.  A  temperature  with  its  accompanying  Maxwell 
distribution  may  thus  be  found  which  has  the  same  impact  ionization  rate  A|  as  that 
for  each  of  the  experimentally  determined  field  distributions.  This  equivalent 
temperature  as  a  function  of  field  strength  is  shown  in  Fig.  8  for  two  different  deep 
lying  acceptors,  Zn  and  Cu.  The  equivalent  temperature  Is  quite  low  and  changes  only 
slowly  with  increasing  field.  It  is  quite  similar  for  both  acceptors,  being  a  little 
higher  for  the  Zn  than  for  the  Cu,  as  expected.  Experimentally  It  was  found^  that  it 
was  somewhat  more  difficult  to  ionize  Cu  than  Zn  impurities.  To  explain  completely 
the  difficulty  of  Ionizing  deep  lying  acceptors  in  Ge,  it  Is  necessary  to  consider 
other  factors  than  Just  the  effective  temperature  of  the  hot  carriers  for  ionization. 
Thus  the  recombination  mechanism  is  an  equally  important  factor.  However,  the  equi¬ 
valent  temperatures  for  ionization  obtained  from  the  above  analysis  are  much  lower 
than  had  been  anticipated  previously  and  must  be  considered  a  major  factor  in  the 
observed  behaviour. 

C.  Hot  Carrier  Mobility 

The  energy  distribution  functions  were  used  to  calculate  the  mobility  of  the  hot 
carriers.  Previously  the  mobility  of  the  hot  carriers  had  been  calculated  using  a 
Maxwellian  distribution  function  with  elevated  carrier  tenperatures  -  a  procedure 
which  had  given  poor  agreement  with  the  experimental  mobilities*.  The  calculations 
made  with  the  new  distribution  function  were  equally  poor,  as  shown  in  Fig,  9.  In 
order  to  explain  why  the  energy  loss  calculations  give  good  agreement  between  theory 
and  experiment  while  the  mobility  calculations  do  not,  the  factors  going  into  each 
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need  to  be  compared.  The  energy  loss  rate  depends  only  on  the  energy  distribution 
function  of  the  hot  carriers.  The  mobility  calculation  depends  on  the  momentum 
distribution  function.  It  Is  the  former  and  not  the  latter  which  is  measured*  and 
the  problem  is  how  to  deduce  the  momentum  distribution  function  f(k)  from  the  measured 
energy  distribution  function  f(t).  f(e)  Is  obtained  from  the  distribution  p(f)  by  a 
density  of  states  factor*  cC  Te*.  The  assumption  was  made  that  the  measured  distribu¬ 
tion  function  was  essentially  isotropic*  and  that  the  perturbation  due  to  the  field 
was  linear*  and  given  by  the  usual  relationship* 


f<£> 


f 

o 


Uhl 
a  «  ' 


(6) 


where  8f/8«  was  evaluated  by  determining  the  slopes  at  various  energies  of  the 
experimental  determined  f(c)  curves.  A  criterion  for  the  validity  of  this  approxima¬ 
tion  is  that  the  perturbation  be  small  compared  to  the  Isotropic  term  f  (which  need 
not  be  the  same  as  the  undisturbed  distribution  at  zero  field  strength).  The 
criterion  can  also  be  approximated  by  the  condition:  v^  ^  ^thermal  ^  When, 
however,  v^r is  calculated  from  the  measured  mobility*  and  from  the 

average  energy  of  the  experimental  distribution  functions*  it  turns  out  that  the  ratio 
Is  of  the  order  of  1/2*  clearly  violating  the  assumption  of  linearity  of  the  perturba¬ 
tion  of  the  distribution  function. 

A  second  difference  between  the  power  loss  and  the  mobility  calculations  is  that 
the  former  depends  only  on  the  optical  mode  collision  time*  while  the  latter  depends 
on  acoustical*  optical*  and  ionized  inpurity  relaxation  times.  However*  we  do  not 
consider  this  difference  to  be  responsible  for  the  difficulties. 

It  Is  of  interest  to  consider  the  reason  for  the  large  ratios  of  v^r^t/  ~therroaj* 
and  the  consequences  of  it.  The  large  values  of  the  ratio  can  be  attributed  to  the 
strong  inelastic  scattering  tn  p-type  germanium*  which  is  just  due  to  the  optical 
phonon  interaction.  Because  of  the  large  energy  loss  rate*  the  carriers  can  not 
become  very  hot*  so  thermal  rema*ns  relatively  low.  The  drift  velocities  however 
become  high  as  the  field  is  increased.  The  carriers  then  have  a  strong  element  of 
streaming  rather  than  diffusive  motion.  The  implication  is  that  the  mon^ntum  distri¬ 
bution  of  the  carriers  is  strongly  displaced  away  from  £  ■  0  in  the  direction  of  the 
field.  An  Interesting  consequence  of  such  non-i 1  nearly  displaced  distributions  is 
that  the  absorption  of  the  !R  light  should  depend  on  the  direction  of  polarization* 
thus  making  the  Ge  dlchroic  when  the  carriers  are  hot.  We  have  looked  for  and 
observed  such  an  effect  and  are  in  the  process  of  investigating  It  in  more  detail. 
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6.  D.  Brown  and  R.  Bray,  Proceedings  of  the  International  Conf.  on  Semiconductor 

Physics  Prague,  (1960),  (Czech.  Acad.  Sciences,  Prague)  P.  82. 


OSCILLATORY  MAGNETORES I  STANCE  OF  GaSb 
W.  M.  Becker  and  H.  Y.  Fan 

Previous  measurements^  on  n-type  GaSb(Te)  showed  that  at  low  temperatures  the 

transverse  magnetoresistance  is  small  for  low  carrier  concentrations,  and  rises 

18-3 

sharply  at  an  nc  ~  1 .2  x  10  cm  where  conduction  in  a  <  111  >-valley  band  sets  In. 
These  measurements  were  carried  out  to  14  kllogauss.  Measurements  made  during  this 
period  extended  the  range  of  investigation  to  30  kllogauss. 

For  samples  with  Fermi  levels  close  to  the  upper  band  edge,  the  measurements 
reveal  single^per lod  Shubnikov-de  Haas  oscillations  in  both  the  transverse  and  longi¬ 
tudinal  magnetoresistance.  The  transverse  and  longitudinal  oscillations  are  In  phase 

for  each  sample,  and  the  period  is  given  by  the  zero  field  Fermi  energy.  Curves 
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showing  behavior  typical  of  samples  having  concentrations  nc  ~  1 .2  x  10  cm  are 
presented  in  Figure  10.  The  an^lltude  of  the  oscillations  increases  sharply  near 

this  concentration.  The  average  component  of  the  transverse  oscillations  increases 
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proportional  to  H  ;  (Ap/pH  )  transverse  also  increases  sharply  it  n(  ~  1.2  x  10  cm  . 

The  results  for  a  number  of  samples  are  summarized  In  Figure  11.  Using  an  effective 
mass  of  0.052  mQ  for  the  center  band,  the  oscillation  at  30  kllogauss  Is  found  to 
correspond  to  quantum  level  ~  10  for  the  samples  Investigated. 

Since  both  the  transverse  magnetoresistance  and  the  oscillation  amplitudes 
increase  sharply  over  the  same  range  of  values  of  the  Hall  coefficient,  the  oscilla¬ 
tions  appear  to  be  correlated  with  the  onset  of  conduction  in  the  higher  lying 
<  ill  >  valleys.  Other  measurements^  show,  however,  that  the  4.2°K  Hall  mobility 
rises  with  concentration,  suggesting  that  the  quantum  condition,  •  t~  1,  may  not  be 
satisfied  at  low  concentrations  for  carriers  in  the  center  band. 
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For  uniaxial  compress  tonal  stress  along  the  <  111  >  direction,  It  Is  expected 
that  there  will  be  a  relative  shift  of  one  of  the  <  111  >  valleys  toward  the  central 
minimum.  If  the  effects  discussed  above  depend  on  the  position  of  the  Fermi  level 
with  respect  to  the  higher  lying  <  111  >  valleys,  the  following  behavior  would  be 
expected:  (1)  a  shift  in  the  oscillations  toward  lower  fields,  (2)  an  increase  In 
the  average  conponent  of  the  transverse  oscillations,  and  (3)  an  increase  In  the 
amplitude  of  the  oscillations.  The  first  two  effects  are  clearly  Indicated  for  a 

3 

sample  with  R^  ^ ° k  "  “  5.49  cm  /coulomb  as  shown  In  Figure  12,  but  no  appreciable 
change  is  noted  In  the  oscillation  amplitude.  The  shift  of  the  oscillations  corres¬ 
ponds  to  a  relative  decrease  In  energy  between  the  center  band  minimum  and  one  of  the 
<  111  >  valleys  of  ^  1  x  10  ^  ev  for  X  ■  3  x  10®  dynes/cm^. 

The  results  on  both  the  transverse  magnetoresistance,  the  oscillatory  behavior, 
and  the  pressure  effects  are  all  consistent  with  a  picture  of  alternating  shift  of 
carriers  from  one  band  to  another  with  changing  field.  However,  calculation  gives 
only  ~  1  shift  of  carriers  from  one  band  to  another  for  the  quantum  level  Involved 
whereas  the  experimentally  observed  variation  is  much  greater  (  x  10  at  30  kilogauss). 
Other  mechanisms  must  therefore  also  be  considered  to  account  for  the  effects. 


1.  W.  M.  Becker,  A.  K.  Ramdas,  and  H.  Y.  Fan,  J.  Appl.  Phys.  32.,  2094  (1961). 


GALVANOMAGNETIC  EFFECTS  IN  HEAVILY  DOPED  P-TYPE  GERMANIUM 

G.  Sadaslv 

The  Hall  effect  and  magnetoresistance  In  several  samples  of  germanium  with 

17  18  “3 

gallium  inpurity  concentrations  ranging  from  6  x  10  to  6  x  10  cm  have  been 
measured.  The  measurements  were  made  at  78°K,  4.2°K,  and  1.3°K;  with  magnetic  fields 
ranging  up  to  30  kfi. 

We  first  discuss  the  results  obtained  at  78°K.  The  values  of  the  resistivity, 
and  the  Hall  coefficient  at  very  low  fields,  for  the  different  samples  are  given  In 
Table  I.  The  direction  of  the  current,  and  the  plane  in  which  the  magnetic  field  was 
rotated  in  making  the  magnetoresistance  measurements,  are  also  Indicated  in  the  Table. 
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TABLE  I  T  -  78°K 


Sample 

R 

3  °  -1 

cm  cou 1 

po 

ohm  cm 

<  i  > 

<  H  > 

Ito 

1  a 

8.46 

.0139 

<  100  > 

<  100  > 

b 

<  110  > 

<  001  > 

II  a 

5.87 

.0106 

<  100  > 

<011  > 

b 

<  110  > 

<  001  > 

III  a 

4.13 

.0081 

<  100  > 

<  100  > 

b 

<  110  > 

<  001  > 

IV  a 

.95 

.0026 

<  100  > 

<  01 1  > 

b 

<  110  > 

<  001  > 

No  longitudinal 

magnetores I  stance 

was  detected  in  any  of  the  sanples.  The 

transverse  magnetoresistance  data  for 

the  sanples  (a) 

are  given 

in  Fig.  13.  Here 

A  p/p  IT  is  plotted  as 

a  function  of  H. 

The  data  for 

the  samples 

(b)  were,  within 

experimental  error,  the  same  as  those  for  the  samples  (a).  As  can  be  seen  A  p/p  is 

2 

approximately  proportional  to  H  ,  though  saturation  effects  are  clearly  apparent  in 
the  purer  samples. 

These  experimental  results  have  been  analyzed  on  the  basis  of  the  following 
model.  The  structure  of  the  valence  band  in  pure  germanium  is  known  from  cyclotron 
resonance  experiments.  We  assume  that  in  the  heavily  doped  material  the  band  struc¬ 
ture  Is  unaltered,  and  that  there  are  no  localized  states,  so  that  we  have  holes 
in  the  normal  valence  band.  There  are  two  types  of  holes,  the  heavy  holes  with 
wf*  ■  .36  mQ,  and  light  holes  with  mi*  *  .044  m^.  The  ratio  of  the  number  of  heavy 
holes  to  the  number  of  light  holes  is  p^/p^  m  24.  The  orientation  dependence  of  the 
magnetoresistance  suggests  that  the  warping  of  the  bands  may  be  neglected,  so  that 
they  can  be  treated  as  simple  spherical  bands.  It  is  also  assumed  that  the  scattering 
mechanism  is  the  same  for  both  types  of  holes,  and  that  the  relaxation  time  t  oc  ts, 
where  c  is  the  energy  of  the  carrier;  the  constant  of  proportionality  is  different  for 
the  two  types  of  holes.  S tuple  transport  theory  can  then  be  used  to  calculate  the 
magnetoresistance.  For  the  sanples  that  we  are  considering  neither  the  non-degenerate 
nor  the  completely  degenerate  approximations  can  be  made  in  calculating  the  transport 
coefficients,  the  exact  Fermi  integrals  were  therefore  used  in  the  computations. 

We  denote  the  ratio  of  the  conductivity  mobilities  for  the  two  types  of  holes  by 
a.  The  best  fit  to  the  data  is  obtained  for  a  value  of  a«3,  and  s  *  3/2.  This  Is 
the  result  one  would  expect  if  the  scattering  were  due  to  ionized  Inpurltles.  The 
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observed  variation  of  Hall  coefficient  with  magnetic  field  is  also  consistent  with 

2 

the  simple  model.  The  values  of  R  as  a  function  of  H  for  sample  II  are  given  In 
Fig.  14. 

The  assumptions  made  Indicate  that  R  pe  m  2.  Actual  chemical  analysis  on  similar 
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specimens  by  Bernard  et.  al.  indicate  that  this  Is  indeed  the  case  If  <10  cm  . 

The  model  is,  of  course,  an  oversimplified  one;  and  the  experimental  results  are 
not  accurate  enough  to  decide  whether  this  Is  the  only  possible  interpretation  of  the 
data.  Measurements  were  therefore  made  at  4.2°K  and  1.3°K.  In  this  case  one  would 
expect  that  the  samples  would  show  complete  degeneracy,  and  a  smaller  magnetoresis¬ 
tance  effect  would  be  seen.  The  experimental  results,  however,  were  completely  at 
variance  with  expectation.  The  transverse  magnetoresistance  at  4.2°K  was  larger  than 
at  78°K;  and  showed  a  linear  dependence  on  field  rather  than  a  quadratic  one.  A 
longitudinal  magnetoresistance  was  also  observed  which  was  comparable  to  the  trans¬ 
verse  one  for  fields  up  to  6  kG,  and  then  showed  signs  of  saturation  at  higher  fields. 
The  results  are  shown  In  Fig.  15.  The  behavior  of  the  samples  (b)  was  the  same  as 
that  of  the  samples  (a). 

The  Hall  coefficient  also  varied  more  rapidly  with  field  at  4.2°K  than  at  78°K. 

2 

The  values  of  R  as  a  function  of  H  for  sample  II  are  shown  In  Fig,  14. 

At  1.3°K  the  magnitude  of  the  magnetoresistance  effect  Increased.  The  results 

for  sanple  II  are  given  in  Fig.  16. 

2 

Keesom  and  Seidel  have  measured  the  heat  capacity  of  heavily  doped  p-type 
material  In  the  temperature  range  from  0.5°K  to  4.2°K.  Their  results  are  in  agreement 
with  the  predictions  of  the  simple  degenerate  gas  model.  Thus  we  must  conclude  that 
though  the  magnitudes  of  Rq  and  pQ  change  only  slightly  in  going  from  78°K  to  4.2°K, 
the  scattering  mechanism  is  different  at  the  two  temperatures. 


1.  W.  Bernard,  H.  Roth,  and  W.  D.  Straub,  Bull.  Am.  Phys.  Soc.  8,  224  (1963). 

2.  P.  H.  Keesom  and  G.  Seidel,  Phys.  Rev.  113.  33  (1959), 


PIEZORESISTANCE  OF  n-GaAs 
R.  J.  Sladek 

Piezoresistance  measurements  were  begun  previously^  on  oriented  single  crystal 

specimens  of  n-type  GaAs  in  order  to  determine  whether  or  not  they  would  indicate 

appreciable  conduction  In  subsidiary  minima  not  located  at  &  m  0  since,  recently, 

2 

magnetoresistance  measurements  of  Kravchenko  and  Fan  revealed  anisotropy  in  magneto- 
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resistance  which  could  be  explained  as  due  to  there  being  appreciable  conduction  In 
(100)  type  minima. 

We  have  continued  the  plezores (stance  investigation  by  measuring  the  piezore¬ 
sistance  of  a  number  of  samples  as  a  function  of  temperature  between  77°K  and  300°K, 

14 

The  samples  had  apparent  carrier  concentrations  at  room  temperature  between  8  x  10 
—3  16  -3 

cm  and  3,7  x  10  cm  and  were  oriented  with  the  long  dimension  in  the  [100],  [110], 
or  [211]  direction.  Uniaxial  compression  and  measuring  current  were  applied  along 
this  dimension  and  the  piezoresistance  was  measured  along  it  except  as  indicated. 

The  longitudinal  piezoresistance  as  a  function  of  temperature  for  a  number  of 
samples  is  plotted  in  Fig.  17  in  which  x  is  actually  the  relative  change  of  resistance 
per  unit  uniaxial  compressive  stress.  (The  numbers  on  the  curves  indicate  the  carrier 
concentration  at  room  temperature,  e.g.  1.5  -  10^  means  1.5  x  10^  carriers/cm3.) 

From  Fig.  17  It  can  be  seen  that:  (1)  x  decreases  with  increasing  carrier  concentra¬ 
tion,  and  (2)  the  samples  fall  into  two  groups  each  having  a  different  kind  of  temper¬ 
ature  dependence  for  x.  One  group  exhibits  a  monotonic  increase  with  decreasing 
temperature  and  the  other  group  exhibits  a  minimum  below  room  tenperature.  The 
position  of  the  minimum  depends  on  carrier  concentration  and  moves  to  higher  tempera¬ 
tures  as  the  carrier  concentration  increases.  If  this  pattern  holds  for  even  higher 
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concentrations  the  negligible  x  measured  by  Sagar  on  a  sample  with  a  somewhat 
greater  concentration  than  our  highest  concentration  is  quite  consistent  with  our 
results. 

The  kind  of  temperature  dependence  exhibited  by  x  can  be  correlated  with  the 
behavior  of  the  Hall  coefficient,  R^,  as  a  function  of  temperature  which  is  plotted 
in  Fig,  18.  (A  field  of  4500  gauss  was  used  for  the  data.)  From  Fig.  18  it  can  be 
seen  that  the  samples  again  fall  into  two  groups.  For  one  group  |  R^  |  increases  very 
greatly  with  decreasing  temperature  and  it  is  these  sanples  which  have  a  large  x 
which  is  monotonic  with  temperature.  For  the  other  group  |  R^  |  increases  only 
slightly  as  temperature  goes  down  and  these  are  the  sanples  in  which  x  goes  through  a 
minimum  as  a  function  of  temperature. 

In  order  to  elucidate  the  origin  of  the  piezoresistance,  the  effect  of  crystal 

orientation  was  studied,  x  was  found  to  be  quite  insensitive  to  crystal  orientation. 

For  example,  for  a  [100]  length  sample  and  a  [110]  length  sanple  each  having  about 
15  3 

1 .8  x  10  carriers/cm  (at  room  temperature)  it  was  found  from  the  longitudinal 
effects  that  x^  -  65  x  10"12  cm2/d  and  ^  (x^  +  xJ2  +  x^)  ■  67  x  10*12  cm2/d  and  that 
from  the  transverse  voltage  with  longitudinal  current  and  stress  mb  x  10"'2  cm2/d. 
Combining  these  results  “  *12  "  1  x  I  O’  cm  /d.  Thus  both  shear  coefficients 
and  X|j  -  x]2  are  very  small  compared  to  .  This  Indicates  that  the  large 
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longitudinal  piezoresistance  effects  are  not  due  to  transfer  of  electrons  between 
multivalley  minima  and  hence  conduction  must  be  overwhelming  In  the  minimum  at  ■  0. 

(Orientation  effect  data  were  also  taken  on  samples  of  the  group  showing  a 
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minimum  In  *  vs  T  with  the  same  type  of  results  as  for  the  1 .8  x  10  cm  samples.) 

Finally  we  measured  the  effect  of  stress  on  the  Hall  coefficient  (at  4500  gauss) 
for  samples  In  both  groups.  In  all  cases 

A  |  Ru  | 

|  RH  |  X  *  *Hal  1  **  *long. 

Thus  the  piezoresistance  seems  to  be  due  to  a  decrease  in  carrier  concentration  with 
increasing  uniaxial  compression.  Probably  the  most  likely  explanation  Is  a  change  in 
the  ionization  of  deep  donor  levels  due  to  the  conduction  band  edge  moving  away  from 
them  due  to  the  applied  stress.  The  two  types  of  temperature  dependence  for  the 
piezoresistance  may  be  a  consequence  of  different  kinds  of  compensation  plus  different 
proportions  of  shallow  to  deep  donor  levels. 

Thermoelectric  power  measurements  are  planned  to  check  this  last  notion  by 
locating  the  Fermi  level.  Piezoresistance  measurements  using  hydrostatic  pressure 
are  also  planned  to  obviate  the  transverse  measurement  needed  when  only  uniaxial 
stress  is  used. 


1.  r.  j.  siadek,  SginKwdiKtw.  RraarsMIghth.  Quarttrly  and  LLaal  Hwrcrt.  Oct.  1, 

1961  to  Sept.  30,  1962. 

2.  A.  F.  Kravchenko  and  H.  Y.  Fan,  Reot.  of  Inti.  Conf.  on  Physics  of  Semlcond. 
Exeter  1962.  p.  737. 

3.  A.  Sager,  Phys.  Rev.  JJ2,  1533  (1958). 

OPTICAL  PROPERTIES  OF  lll-V  SEMICONDUCTOR  COMPOUNDS 
E.  J.  Johnson  and  H.  Y.  Fan 

Infrared  measurements  at  1 .5°K  show  a  step  on  the  absorption  edge  of  InSb. 
Under  applied  magnetic  field  the  step  develops  Into  a  peak  which  moves  to  higher 
photon  energy  and  splits  into  two  consonants*.  This  behavior  can  be  described 
approximately  by: 


(hv)+  w  (Eg  -  «)  +  £  (fi  9C  ±  gc  H) 
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where: 

1)  (hv+f  (hv)_  are  the  photon  energies  of  each  peak, 

2)  €  -  .006  ev 

3)  odc  is  the  cyclotron  frequency  of  electrons. 

4)  g  is  the  g-f actor  for  spin  spl \ tt Ing  in  the  conduction  band. 

c  1 

In  analogy  to  previous  work  on  GaSb  the  absorption  was  attributed  to  transitions 
involving  an  Ionized  acceptor  level  and  the  conduction  band.  Related  bound  exciton 
absorption  may  be  present,  particularly  with  applied  magnetic  field.  Recent  measure¬ 
ments  (Fig.  19)  demonstrate  that  the  strength  of  the  absorption  is,  indeed,  related 
to  impurity  concentration.  Further,  the  following  considerations  show  that  the 
impurity  level  involved  is  one  near  the  valence  band.  1)  The  threshold  for  the  absorp¬ 
tion  occurs  at  hv  <  Eg  -  .006  ev.  The  related  ionization  energy  Is  then  |  Ej  |  >  .006 

ev.  Therefore  the  level  lies  between  E^  -  .006  and  EyB  +  .006.  2)  The  absorption 

is  observed  in  degenerate  n-type  material,  therefore,  with  the  Fermi  level  in  the 
conduction  band.  3)  The  absorption  is  observed  in  p-type  samples  with  the  Fermi  level 

near  E^^  +  .006  ev.  Any  level  between  E^  and  E^^  +  .006  ev  will  be  In  different 

charge  states  in  the  respective  samples  and  cannot  give  rise  to  the  same  absorption. 
Therefore  the  level  Is  near  the  valence  band  -  either  an  ionized  acceptor  or  a  neutral 
donor.  The  level  is  probably  an  ionized  acceptor,  possible  Zn. 

The  absorption  edges  of  various  other  lll-V  compounds  and  silicon  have  been 
surveyed  at  1 .5°K  for  bound  exciton  and  impurity  absorption. 

GaAs  shows  a  step  in  absorption  at  about  1.495  ev  which  is  probably  related  to 
an  impurity  level  of  about  .023  ev.  This  step  is  not  significantly  affected  by  mag 

2  3 

netic  fields  as  high  as  22  kG,  This  step  in  absorption  has  been  observed  by  others 

3 

Nathan  sees,  in  addition,  a  peak  In  emission  apparently  associated  with  the  reverse 
transition. 

InP  has  a  slowly  rising  tall  on  the  absorption  edge  of  about  the  same  magnitude 
at  both  room  tenperature  and  1 ,5°K.  At  the  lower  temperature  several  steps  occur 
which  are  insensitive  to  a  magnetic  field  of  22  kG.  The  steps  are  probably  associated 
with  Inpurlty  absorption.  The  absorption  rises  sharply  to  beyond  10,000  cm  *.  Very 
thin  samples  are  necessary  for  measurement  at  shorter  wavelengths. 

InAs  displays  a  sharply  rising  edge  from  an  absorption  of  2  cm" *  to  about  1000 
cm"*.  This  may  be  the  Intrinsic  absorption  edge  or,  more  likely,  the  threshold  for 
impurity  or  bound  exciton  absorption. 

The  absorption  edge  In  silicon  begins  with  Indirect  transitions  and  should 
present  an  excellent  situation  to  study  bound  exciton  transitions.  These  can  take 
place  without  phonon  participation  and  appear,  therefore,  In  a  region  of  the  spectrum 
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where  no  band  gap  transitions  occur.  Haynes  has  observed  strong  emission  peaks  In 
silicon  which  he  attributes  to  excltons  bound  to  neutral  Impurities.  Attempts  were 
made  to  observe  the  corresponding  absorption.  Samples  with  Impurity  concentrations 
(B)  as  high  as  3.3  x  10  cm  were  examed.  Ho  such  absorption  was  observed  although 
the  background  absorption  was  less  than  .02  cm  \  Samples  doped  with  B  and  Sb  com¬ 
pensating  Impurities  (Np  «  1,5  x  10^7,  «  8  x  10^)  have  also  been  studied.  Again, 

It  has  not  been  possible  to  observe  absorption  due  to  excitons  bound  to  ionized 
impurities. 

l#  E.  J.  Johnson  and  H.  Y.  Fan,  Semiconductor  Research,  Signal  Corps  Contract 
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INFRARED  ABSORPTION  OF  OXYGEN  IN  SILICON  AND  GERMANIUM 
UNDER  HIGH  RESOLUTION 
A.  K.  Ramdas  and  R.  L.  Aggarwal 

Oxygen  enters  Into  silicon  crystals  grown  from  silicon  melted  in  quartz  crucibles. 

An  absorption  band  in  such  samples  at  around  9  microns  has  been  correlated  to  dispersed 

l  2 

oxygen  in  silicon  by  Kaiser,  Keck,  and  Lange  .  Hrostowski  and  Kaiser  established 
that  it  is  due  to  the  ant i-symmetr ic  mode  of  vibration  of  the  bent  Si-O-Si  unit.  An 
analogous  band  in  germanium  is  observed^  in  germanium  at  11.7  microns.  At  temperatures 
below  ~  150°K  the  9  micron  band  in  silicon  exhibits  a  remarkable  fine  structure^.  The 
object  of  the  present  investigation  is  to  examine  the  fine  structure  for  the  above 
mentioned  bands  under  high  resolution  as  there  is  considerable  Interest  regarding 

3 

their  origin  . 

Crucible  grown  samples  of  silicon  with  no  intentional  doping  have  been  used. 

Oxygen  was  intentionally  introduced  in  germanium  via  the  ambient.  A  high  resolution 
grating  instrument  equipped  with  blazed  gratings  was  employed.  The  effective  resolu¬ 
tion  under  which  the  samples  were  measured  was  ~0.5  cm  The  absorption  spectra  were 
measured  using  a  glass  optical  cryostat  and  liquid  helium  as  coolant.  The  samples  were 
cooled  by  conduction  and  the  sample  temperature  is  estimated  to  be  10-15°K. 
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Figure  23  shows  the  absorption  spectra  for  silicon.  Four  absorption  peaks  have 
been  clearly  observed  at  1128.1g,  ll32.7Qt  1134.5^  and  1136.3^  cm  '  respectively. 

Weak  peaks  have  also  been  observed  at  1 1 24 . 7^ ,  1126.5p  and  1 1 30 . 9Q  cm  '  respectively. 
The  1 1 28 •  1 0  cm"'  peak  shows  a  variation  in  strength  from  sample  to  sanple.  A  number 
of  other  sharp  but  weak  peaks  have  been  observed  In  the  range  1030*1120  cm  at 

1084.92,  1099 .3^,  1103. 9^  and  1107.6^  cm"'.  The  origin  of  these  bands  is  not  yet 
clear. 

3 

According  to  Hrostowski  and  Alder  oxygen  bonding  the  two  nearest  neighbour 
silicon  atoms  in  the  SI-0-S1  unit,  can  tunnel  from  a  potential  minimum  to  any  one  of 
the  other  symmetrically  equivalent  minima,  the  symmetry  resulting  from  the  fact  that 
the  silicon  atoms  bonded  to  the  oxygen  atom  are  along  <  111  >axls.  The  harnonlc 
oscillator  levels  belonging  to  the  ant i-symnetr ic  mode  split  as  a  result  of  the 
tunneling.  The  observed  fine  structure  and  its  temperature  variation  can  be  qualita¬ 
tively  understood  on  this  model.  In  the  present  work  the  fine  structure  has  been 
observed  with  all  its  details  clearly  resolved  and  hence  ft  will  be  possible  to 
follow  the  temperature  variation  of  the  intensity  of  the  components  and  deduce  a  model 
for  the  tunneling  phenomenon.  An  optical  cryostat  with  an  exchange  gas  system  has 
been  set  up  for  this  purpose. 

The  ant i -symmetric  mode  of  vibration  of  Ge-0-Ge  unit  in  germanium  responsible  for 
the  11.7  micron  band  also  exhibits  fine  structure  at  low  temperature  and  under  high 
resolution  (See  Figure  24).  The  strongest  peak  at  862.20  cm"'  is  followed  by  peaks 
at  862. 8^  and  863.8  cm  '  respectively.  It  is  interesting  to  note  that  unlike  in  the 
case  of  silicon,  the  low  frequency  conponent  is  the  strongest.  A  detailed  examination 
of  the  fine  structure  Is  also  planned  as  its  origin  appears  to  be  similar  to  the 

analogous  case  of  silicon  discussed  above. 

4 

It  is  known  that  oxygen  in  silicon  on  heat  treatment  in  the  temperature  range 
300-600°C  forms  donor  complexes.  A  series  of  absorption  bands  have  been  observed  in 
a  heat  treated  sanple  at  1043. 8^,  1052. 0^,  1057. 6Q,  1059. 7^,  and  1090. 1Q  cm"'  and 
are  presumably  due  to  vibrations  of  various  silicon  oxygen  conplexes  which  behave 
like  donors  (See  Figure  23).  Some  of  these  have  been  reported  by  Kaiser,  Hrostowski, 
and  Thurmond^.  The  donor  conplexes  show  electronic  excitation  spectra  at  low  temper¬ 
atures^9^,  and  it  appears  that  there  are  many  species  of  such  donors,  it  is  of  inter¬ 
est  to  correlate  the  vibrational  absorption  bands  with  the  electronic  excitation  bands. 
Further  work  is  in  progress  along  these  lines* 
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II.  IRRADIATION  EFFECTS 

ELECTRON  PARAMAGNETIC  RESONANCE  IN  SILICON  AND  GERMANIUM 
G.  S.  Newell,  W.  Jung,  D.  Trueblood,  end  E.  Hale 


Introduction 

In  another  report  (hereafter  referred  to  as  "A") *  we  summarized  the  annealing 
curves  and  spin  Hamiltonians  of  the  paramagnetic  defects  which  have  been  Isolated 
during  the  annealing  of  plle-lrradlated  high  resistivity  silicon.  We  here  discuss 
their  hyperfine  structures  and  also  derive  a  better  value  for  the  absolute  production 

rates.  The  material  summarized  in  these  two  reports  Is  discussed  in  detail  in  a 

2  3 

thesis  by  W.  Jung  and  has  been  submitted  for  publication  . 

Helium  tenperature  investigation  of  germanium  irradiated  with  electrons  at  77°K 
Is  in  progress. 

HyMrfl,.ng-?tr.,UGLyrs 

For  Centers  (II,  III),  (V,  VI),  and  N,  which  were  observed  In  large  enough  con~ 
centrations,  the  weak  hyperfine  satellites  arising  from  the  magnetic  interaction  of 
the  electronic  magnetic  moments  with  neighboring  Si  nuclei  (I  w  4. 7f  abundant) 
were  observed.  The  recorder  trace  of  a  Center  (II,  III  spectrum  (H  |  <  1 10  >)  which 
is  reproduced  in  Fig.  20  shows  these  satellites.  For  Center  (II,  III),  their  ampli¬ 
tude  ranges  from  4  to  6^  of  satellites.  For  Center  (II,  III),  their  amplitude  ranges 
from  4  to  6^  of  the  main  lines  and  is  well  above  the  noise  level.  The  separation  from 
the  center  line  varies  from  23  to  37  gauss  and  is  fairly  well  resolved.  For  Centers 
(V,  VI),  and  N,  however,  the  hyperfine  signals  appear  to  be  much  weaker  (only  2  to  3^) 
and  are  overshadowed  by  other  centers  at  most  angles,  while  the  separation  seems  to  be 
comparable  to  that  of  Center  (II,  III). 

In  addition,  some  of  the  well  isolated  branches  of  Centers  (II,  III),  (V,  VI), 
and  (VII,  VIII)  show  much  more  closely  spaced  satellites  ("triplet  structures")  (Fig. 
20)  arising  from  hyperfine  interactions  with  more  distant  nuclei.  The  lines  on  each 
side  are  separated  from  the  center  line  by  about  3  gauss  and  the  amplitudes  are 
20  to  29^  of  the  center  lines. 

The  presence  of  several  centers  in  each  spectrum  and  an  unfavorable  signal  to 
noise  ratio  prevented  determination  of  the  hyperfine  tensor  with  any  reliability 
except  for  Center  (II,  III).  Even  for  Center  (II,  III)  the  "triplet"  lines  were  so 
obscured  by  other  branches  that  the  hyperfine  tensor  for  those  satellites  was  not 
determined. 

4 
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By  Including  the  hyperfine  interaction  term  in  the  spin  Hamiltonian,  the  hyper- 
fine  satellites  are  described  by 

>f-  0  H-q/$  +  S'O-S  +  1 -A-S  +  0  g  H-l.  (1) 

^  **M*m*  ******  ******  *n  *n 


with  $  ■  1 . 

Neglecting  the  last  term,  the  first  order  energy  shift  due  to  one  magnetic 
nucleus  is 


E(m^m|)-E(msf0)  »  n^mj  |jk  •  Aj  «  msm|Ae»  (2 ) 

2  2 

where  A  ■  k  •  A  •  k;  m  ,  m,  are  the  electronic  and  nuclear  magnetic  quantum  numbers; 
e  ~  **  s  i  * 

and  k  Is  the  unit  vector  in  the  direction  of  the  effective  field:  k  •  (g  H)~  g  •  H. 

**  **  c  **  ** 

As  the  microwave  transitions  correspond  to  Am^  ■  1,  Am|  ■  0,  the  separation,  AH, 
between  the  main  and  satellite  lines  is  given  by 

-  m|Ae  (3) 

1  29 

where  m  +  j  f°r  SI  •  Hence  the  hyperfine  interaction  tensor,  is  determined 
from  the  angular  variation  of  the  quantity 

A2  -  (2  g^gAH)2  (4) 

by  applying  the  least  squares  procedure  developed  for  the  fine  structure  tensors. 

The  second  order  corrections  to  the  hyperfine  splittings  are  obtained  from  the 
Hamiltonian  (Eq.  ))  by  a  perturbation  treatment  in  which  T  is  kept  in  operator  form*. 
The  result  is  an  effective  nuclear  Hamiltonian  which  for  I  ■  ^  takes  the  form 

U±  "  +  V|mr'» 


{C3m2  -  $($  +1)]  [ic*o*A  -  (TT‘9-1c)  U*A] 

-  j  +  1)  -  m2]  (Oet  Aj)  lc*A|” +  g^ff 


with 
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Here  9e  *  1^  ‘  $)  *s  the  effective  g  value,  a  Is  the  zero  field  splitting  tensor,  and 
Det  A  Is  the  determinant  of  A.  The  last  term  is  the  direct  interaction  of  the  nucleus 
with  the  external  field:  It  Is  a  first  order  term,  but  has  about  the  same  magnitude 
as  the  second  order  corrections.  The  significant  feature  is  that  is  a  vector 

which  depends  on  |m|,  not  m,  and  hence  aids  the  field  seen  by  the  nucleus  in  one 
|m|  ■  1  level  and  opposes  it  in  the  other.  The  corrected  hyperfine  splittings  become: 


AH 


1 

*eP 


where  7.,  7  are  the  values  of  7»  i  in  the  m  *  +  l,  and  m  ■  0  levels.  For  Center 
l  o  |m|  s  "  s 

(11,111)  at  k-band  the  corrections  are  comparable  to  experimental  error  (~  0.5  gauss) 
and  are  neglected.  Under  more  favorable  conditions,  however,  it  may  be  possible  to 
use  the  difference  in  the  hyperfine  splittings  to  assign  absolute  m$  values  to  the 
spectra. 

The  observed  and  calculated  angular  variation  of  the  hyperfine  splitting  of 
Center  (II,  III)  are  plotted  in  Fig.  21.  To  measure  a  splitting,  both  the  satellite 
and  the  main  line  must  be  visible;  thus  in  constructing  Fig.  21  points  were  taken  from 
any  of  the  14  main  branches  which  happened  to  be  usable  and  are  keyed  to  the  corres¬ 
ponding  branches  plotted  in  Fig.  22  by  using  similarly  coded  lines.  For  Center  (II, 
III),  the  following  values  are  obtained: 


A,  -  (44. 1  +  0.4)  x  Kf4  cnf 1  (6  -  36.5°  ±  0.5°) 

A2  -  (43.6  +  0.4)  x  10“4  cm-1 

A3  -  (68.8  +  0.4)  x  I0-4  cm"’ 

where  A^  lies  along  <  110  >  and  the  axes  are  defined  as  usual^.  These  show  a  nearly 
axial  <  111  >  symmetry;  A|j  ■  68.8  x  10~^  cm  *  and  Aj_«  (Aj  +  A^)  *  43.9  x  10~^  cm"* 
with  A n  lying  1.2°  from  <  111  >  toward  <001  >  in  the  (110)  plane.  In  order  to 
examine  the  overall  fit  of  the  parameters,  the  angular  variations  of  the  hyperfine 
satellite  lines  were  calculated  from  the  fitted  values  of  a,  and  A,  and  plotted 
for  all  of  the  28  hyperfine  branches  corresponding  to  the  14  main  lines.  Figure  22 
shows  the  actual  angular  variation  plot  for  the  unannealed  Sample  2*.  Exactly  the 
same  curves  are  obtained  for  Sample  4®.  The  thick  lines  represent  the  14  main  lines 
and  the  thin  lines  the  corresponding  hyperfine  satellites,  with  both  the  parent  and 
satellite  lines  coded  in  the  same  way  as  in  Figs.  10  to  14  of  (A).  The  fitted  curves 
account  for  all  of  the  visible  portions  of  the  hyperfine  satellites  remarkably  well, 


-  20 


In  fact  within  a  line  width*  Further,  ail  the  lines  visible  at  room  tenperature  are 
accounted  for.  Those  points  in  Fig.  22  which  do  not  fall  on  the  calculated  curves 
belong  either  to  the  next  nearest  neighbor  h  f  s.  ("triplet  structure") ,  or  to 
Center  IX* 

The  ratio  of  the  anplitudes  of  the  hyperflne  satellites  to  the  main  lines  was 
measured  at  those  angles  where  both. lines  were  well  Isolated  from  others.  The  observed 
ratio  of  5,12jf  (o  ■  0,78%  from  32  measurements)  is  in  satisfactory  agreement  with 
the  value 


[21  +  l]”1  [2C,  P(l-P)]  /  r2CQ  (1-P)2]  -  0.0494 
29  1 

expected  for  two  equivalent  St  nuclei  having  I  ■  and  P  •  0.047,  and  is  easily 
distinguishable  from  the  values  of  2A7%  expected  for  only  one  nucleus,  7  A2%  for 
three,  and  9.89  i  for  four. 

For  a  spin  I  center  such  as  Center  (II,  III),  the  orbitals  localized  on  the  two 

atoms  overlap  very  slightly,  giving  rise  to  an  exchange  coupling.  The  observed  hyper- 

29 

fine  satellites  correspond  to  the  case  In  which  only  one  of  the  two  nuclei  Is  St 
(I  ■  ^),  the  other  being  Si2®  (I  *  0).  Within  the  triplet  manifold,  the  appropriate 
hyperflne  Hamiltonian  Is 


H.  .  ■  i  I  -a  -S, 

hfs  2  ~|  ~i  ~1 


+  2  if  ~l 


where 


t* 


2  V 


I  -  1  or  2. 


The  hyperflne  splitting  of  the  energy  level, 

E(ms,m,)  -  E(ms,0)  -  |  k  •  Aj  |  -  £  |  Jk  •  a^  |, 

is  exactly  the  same  as  for  S  •  2,  the  value  of  ms  ■  £  I  canceling  the  factor  of  j. 

The  fact  that  a  single  set  of  satellites  of  5)f  relative  amplitude  is  observed  indi¬ 
cates  that  £]  "  *nd  that  the  two  nuclei  are  accurately  equivalent.  This  suggests 
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that  tf  the  distortion  of  the  hfs  axis  away  from  <  111  >  is  due  to  the  mutual  repul¬ 
sion  of  the  dangling  bonds,  these  bonds  are  anti-parallel  rather  than  parallel;  In 
the  latter  case  the  hfs  axis  would  differ  by  2.4°  and  the  nuclei  would  be  clearly 
inequivalent* 

7  8 

Watkins  and  Corbett  ’  discussed  the  hyperfine  interaction  of  the  electron  irra¬ 
diation  centers  in  silicon  (all  S  ■  ~)  on  the  following  basis: 

The  wave  function  of  the  electron  Is  approximated  by  3s  and  3p  orbitals  localized 
at  various  lattice  sites  J: 


♦  ■ z)  ’b  <y*J+V3pJ)- 


2  A 

where  a.  and  p.  are  normalized  according  to  a.  +  p/ 

J  J  2  J  4 

fractional  contribution  from  site  J,  and  a.  and  p  2 

o  7  J  9  J 


■  1*  The  factor  7)j  gives  the 
give  the  relative  weights  of  3s 
and  3p  orbitals.  The  values  of  a,*,  and  p.‘  can  be  calculated  from  the  observed 


A  tensor: 


■j  *  j 


j 


*j  -  3  <A|/+2A1J>  ■  ^  Wh  W  •♦Ss  <°>l. 

bJ  “  3  '  AJ^  "  S  <  r  3  *3p  • 


where  gN  is  the  nuclear  |  factor,  PQ  the  Bohr  magneton,  and  p^  the  nuclear  magneton. 

Using  the  value  |t~  (0)  |  /  <  r  ^  >w  ■  1.4  estimated  from  the  tabulated  Hartree 

i  ,2  24  -3 

functions  and  the  value,  |  +3S (0)  I  “  24  x  10  cm  ,  obtained  from  previous  work  on 
the  Si -A  center,  the  authors  quoted  obtain  the  values 


-  Gtj^T)j2  x  1040  x  I0-4  cm-1 


and 


bj  "  x  36  x  10-4  cm_l  • 

A  similar  analysis  of  the  A  tensor  of  Center  (11,111)  gives 

a  -  52.15  x  10-4  cm'1,  b  -  4. 15  x  I0-4  cm"1 
and  hence  O2  -  0.176,  p2  -  0.824,  and  Tl2  -  0.281. 
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In  terms  of  this  simple  model,  the  result  is  interpreted  as  follows:  two  orbitals 
of  18^  3s  and  62%  3p  character  localized  on  two  neighboring  atoms  account  for  56^  of 
the  total  wave  function  of  each  electron,  the  enhanced  p-like  character  over  the  normal 

3 

25^  3s  and  75^  3p  (sp  )  tetrahedral  orbital  suggesting  that  the  atom  Is  pulled  away 
from  the  normal  site  by  Its  neighbors.  The  rest  of  the  wave  function  is  spread  over 
more  distant  atoms. 

The  side  lines  of  the  '’triplet'*  structure  of  Center  (llf  III)  may  account  for  a 
large  portion  of  the  remaining  wave  function.  The  relative  amplitude  varies  from  20 
to  29^  (23.6ft,  a  ■  3.7f(  from  30  measurements),  and  it  is  difficult  to  distinguish 
among  the  values  of  19.8^  expected  for  eight  equivalent  nuclei,  24.7ft  for  ten,  and 
29.7ft  for  twelve.  As  the  approximately  ten  neighboring  atoms  are  almost  certainly 
not  equivalent  for  arbitrary  orientations  of  the  magnetic  field,  these  satellites 
should  split  Into  groups  characteristic  of  the  structure  of  the  center.  Unfortunately, 
the  resolution  is  such  that  the  lines  merely  broaden  and  disappear.  The  fact  they 
are  visible  over  quite  a  range  of  angles  reflects  the  small  value  of  the  coupling 
constant  for  3p  states. 

Obscuration  of  the  lines  and  inadequate  resolution  thus  prevent  determi nat ion  of 
the  hyperfine  tensor  for  these  atoms  and  only  rather  broad  limits  can  be  placed  on 
their  contribution  to  the  total  wave  function.  The  separations  from  the  center  lines 
are  about  3.1  gauss  and  are  nearly  constant  over  the  region  where  they  are  visible. 

„  ]  _4 

The  separation  corresponds  to  5.8  x  10  cm  Assuming  A  5.8  x  10  cm  as  a 

3  2 

reasonable  guess,  and  assuming  little  disturbed  normal  sp  orbitals  (a  -  0.25  and 
2  2 

0  ■  0.75),  A 1 1  »  1.35  Aj^  and  hence  7]  ~  0.026.  Thus  by  taking  ten  nuclei,  the  contri¬ 

bution  to  the  total  wave  function  amounts  to  26f{.  However,  the  data  are  insensitive 

2  2 

to  the  p  component.  Taking  a  -  1  as  the  extreme  case,  11  »  .006  and  the  contribution 
to  the  wave  function  would  be  only  6# 

The  interpretation  that  the  triplet  structures  arise  from  the  interaction  with  a 

number  of  distant  nuclei  is  confirmed  by  the  occasional  detection  of  a  second  pair  of 

very  faint  satellites  having  the  separations  and  amplitudes  predicted  for  interactions 
29 

involving  two  Si  nuclei  among  the  8  to  12  neighbors. 

The  spectra  of  both  main  and  hyperfine  satellite  lines,  and  their  relative 
intensities,  remain  unchanged  between  300°K  and  77°K.  This  suggests  that  the  elec¬ 
tronic  configuration  is  rather  firmly  locked  and  not  affected  in  this  temperature 
range. 

In  contrast,  the  spectrum  of  the  H  Center  goes  through  a  transition  at  tempera¬ 
tures  between  77°K  and  166°K.  While  the  hyperfine  tensor  at  300°K  cannot  be  deter¬ 
mined  accurately  because  of  an  unfavorable  signal  to  noise  ratio,  the  hyperfine 
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splitting  Is  only  one-half  that  of  the  77°K  spectra,  and  the  relative  amplitude  of 
the  satellites  remains  unchanged  at  about  2%.  This  suggests  that  at  300°K,  the  elec¬ 
tronic  wave  function  Is  distributed  by  a  motional  effect  over  two  atoms  which  are  not 
quite  equivalent.  The  hyperflne  Interactions  Involving  two  similar  but  not  equivalent 
sites  should  give  rise  to  twice  as  many  satellite  lines  as  for  two  equivalent  sites, 
each  having  2.5^  amplitude  and  comparable  separations,  and  would  make  the  experimental 
determination  of  the  tensors  much  more  difficult.  This  may  very  well  be  the  case  of 
Center  (V,  VI)  also;  In  fact,  the  anisotropy  of  the  g  tensor  of  the  center  is  consis¬ 
tent  with  a  pair  of  two  <111  >  dangling  bond  orbitals  oriented  toward  the  center  of 
a  tetrahedron  and  thus  supports  the  suggestion. 

In  summary,  the  h.f*s.  data  unequivocally  confirm  the  model  for  Center  (II,  III) 
consisting  of  two  electrons  localized  In  parallel  (or,  more  likely,  antiparallel) 

3 

dangling  sp  orbitals  centered  on  different  atoms.  It  Is  tenpting  to  suppose  that 
Center  (II,  III)  represents  the  neutral  state  of  the  divacancy,  as  Is  suggested  by  the 
number  of  next  nearest  neighbors,  but  this  assignment  is  incompatible  with  the  work 
of  Watkins  and  Corbett  on  the  Si-C  and  SI-J  centers  on  the  ground  that  Center  (II,  III) 
anneals  at  too  low  a  temperature  and  fails  to  show  re-orientation  near  77°K. 

The  higher  temperature  S  «  1  centers  are  probably  similar  to  Center  (II,  III) 
but  with  the  addition  of  oxygen  bridges. 

Matl  £9ns9ntr?tl9n?: 

Figures  8  and  9  of  (A)  plot  the  amplitude  of  the  fine  structure  lines  of  the 

various  centers  relative  to  100  for  the  Internal  ruby  standard.  Our  original  derlva- 
g 

tion  of  absolute  defect  concentrations  from  these  curves  of  necessity  omitted  cor¬ 
rections  for  the  hyperflne  intensities,  and  in  addition,  contained  an  error  of  a 
factor  of  4  in  the  transition  matrix  element  for  ruby.  A  corrected  treatment  follows: 

From  the  amplitude  of  a  center  relative  to  that  of  the  Internal  ruby  standard, 

Its  absolute  concentration  can  be  estimated  using  the  relation 


n 


s 


T)  a  C  A  (AH  r 
N  r  r  s  s  s 

r  T\  a  C  A  (AH  )2 
s  s  r  r  r 
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V 
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-  5.25  x  10 
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where  the  symbols  are  defined  as  follows: 
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rf 2  rW,Y  *  the  relative  filling  factor  ■  0.216, 
Hrf  sample 

S(S+1  )-m  (m  +1 ) 

transition  probability  factor  ■  - -  *  ■  . 

2S  +  1 


l  l  +3 

-  1,02  instead  of  1  for  the  ms  "  "  2  to  +  2  trans*t*on  Cr  in  ruby  due  to 

the  matrix  element  correction. 


AH^  *  1  ine  width 

AH  ■  line  width  of  ruby  signal  *  15.9  gauss.  (Full  width  between  inflections), 

r  +3  1 6 

■  number  of  Or  spins  In  1.79  mg  of  ruby  m  2.07  x  10  , 

■  sanple  volume  ■  1,03  x  .43  x  .093  *  4.12  x  10~*  cm\ 

A  /A  ■  amplitude  ratio 

s  r 

Cs,  Cp  ■  line  shape  factors:  ■  C^,  assuming  the  same  line  shape  (gauss l an). 

The  factor  of  12  results  from  the  fact  that  the  amplitudes  plotted  in  Fig.  5  are  for 
a  single  branch  normalized  to  unit  multiplicity. 

The  peak  concentrations  of  the  centers  at  various  stages  of  annealing  are  col¬ 
lected  in  Table  II.  The  table  includes  corrections  for  (1)  hyperfine  satellite 
lines  and  for  (2)  singlet  level  populations  for  the  spin  1  centers.  For  Center  (llf 
III),  for  example,  weak  satellites  (4.9^0  and  the  "triplet**  satellites  (23.6$)  give 
rise  to  ~  57*  hfs  correction.  For  Centers  (V,  VI)  and  (VII,  VIII),  only  the  "triplet" 
satellites  are  taken  into  account  and  hence  the  50jt  correct  ion  is  only  approximate. 

For  the  N  center,  no  "triplet"  structure  was  observed  and  the  correction  represents 
the  contribution  of  the  weak  satellites  alone.  The  singlet  level  population  correc¬ 
tions  for  spin  1  centers  were  assumed  to  be  +  |  by  taking  a  statistical  weight  of  1/4 
for  the  singlet  levels  and  neglecting  the  Boltzmann  factor. 

The  value  of  2.07  x  10^  spins  in  the  ruby  standard  is  calculated  from  the  nominal 
concentration  of  O.l^tCr  .  Although  it  is  certainly  in  error,  it  is  being  used  until 
we  are  better  satisfied  with  the  accuracy  of  our  own  determinations.  These  have  been: 

(1)  comparison  of  the  resonance  anplitude  with  that  of  a  known  sample  of  Cu  S0A  • 

1 6  ** 

5  HgO,  giving  I .2  x  10  spins  for  the  ruby,  and  (2),  a  measurement  of  the  optical 

density  at  selected  wavelengths^,  giving  1.4  x  10*®  spins: 

Erratum: 

In  (A),  p.  15,  the  equation  relating  the  observed  fields  for  resonance  to  the 
values  of  g  and  o  read  from  the  calculated  curves  should  read. 


TABLE  II 
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H 


Taking  the  average  elastic  scattering  cross  section  of  silicon  atoms  for  fast  neutrons 

1 7  3 

to  be  «  3  barns,  the  number  of  collisions  Is  2.1  x  10  per  cm  for  a  flux  of  1,4  x 
18 

10  nvt,  giving  peak  production  rates  of  0.01  to  0.05  centers  per  collision.  In  view 
of  the  large  number  of  knock-ons  U  20)  expected,  this  small  value  appears  significant. 


1.  Semiconductor  Research,  Contract  DA  36-039-sc-87394,  Eighth  Quarterly  and  Final 
Report,  October  1,  1961  -  Sept.  30,  1962. 

2.  W.  Jung,  Thesis,  Purdue  University,  1963. 

3.  W.  Jung  and  G.  S.  Newell,  Phys.  Rev.  (to  be  published). 

4.  The  derivation  is  somewhat  lengthy.  (See  reference  (1). 

5.  A2  is  along  <  110  >,  A^  lies  in  the  plane  containing  <001  >  and  <  TlO  >,  making 

an  angle  0  with  <  001  >,  and  A^  lies  in  the  same  plane,  perpendicular  to  Ap 

6.  Sample  2  «  crucible  grown;  Sample  4  ■  floating  zone.  See  table  II  of  (A)  for 
details. 

7.  G.  D.  Watkins  and  J.  W.  Corbett,  Phys.  Rev.  JH,  1001  (I960. 

8.  G.  D.  Watkins  and  J.  W.  Corbett,  Discussions  Faraday  Soc.  3J_,  89  (1961). 

9.  Semiconductor  Research;  Contract  DA  36-039-sc~87394,  Fifth  Quarterly  Report, 
October  1,  1961  -  December  31,  1961,  p.  12. 

10.  The  absolute  absorption  coefficient  of  ruby  vs.  wavelength  is  given  by  T.  H. 
Maiman  et  al ,  Phys.  Rev.  J23,  1151  (1963). 

11.  D.  J.  Hughes  and  J.  A.  Harvey,  Neutron  Cross  Sections.  BNl  325,  Brookhaven 
National  Laboratory,  1955,  The  total  cross  section,  Oy,  varies  from  2  to  4  barns 
for  0.1  to  10  MeV  neutrons. 
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III.  LOW  TEMPERATURE  STUD  I ES 

SPECIFIC  HEAT  OF  LEAD  BETWEEN  0.3  AND  4°K 
P.  H.  Keesom  and  B.  J.  C.  van  der  Hoeven,  Jr. 


It  Is  found  that  C  ,  the  electronic  contribution  to  the  specific  heat  of  super- 

1 

conductive  lead,  below  IK,  exceeds  the  value  calculated  by  the  BCS  theory  by  several 

3 

orders  of  magnitude.  Instead,  C  shows  reasonable  agreement  with  a  T  dependence  for 

cs 

values  of  t  ■  T  /T  between  2  and  22. 
c 

A  lead  sample  of  1  .293  moles  (99 .9995C  pur  I ty)  was  measured  and  the  results  In 
the  normal  and  superconduct Ive  states  can  be  represented  by 


Cn  -  2.98  T  +  1.60  T3  +  0.028  T5  mJ/mole  deg  ,  (I) 

Cs  -  0.144  T3  +  1.60  T3  +  0.028  T5  mJ/mole  deg  .  (2) 

The  accuracy  above  1°K  Is  1  -  2%,  and  2-3^  below  this  temperature. 

The  last  two  terms  In  both  equations  originate  from  the  lattice  vibrations  and 

give  for  the  Debye  parameter  at  0°K:  $  -  106.7  +  0.5°K.  This  Is  in  good  agreement 

o  2  *" 

with  results  calculated  from  velocity  of  sound  measurements  (105  K)  .  The  linear 
term  in  eq.  (1)  arises  from  the  electronic  contribution.  This  agrees  with  the  value 
for  the  coefficient  of  3.06  ±  0.04,  deduced  from  critical  field  measurements  by  Decker, 

3 

et  al.  ,  within  the  combined  limits  of  accuracy. 

The  electronic  contribution  in  the  superconductive  state,  C  ,  is  calculated  in 

es 

the  BCS  theory  to  be: 


C  H  T  -  8.5  exp  (-  1 .44  T  /T)  .  (3) 

6S  C  C 

Decker  et  al.  already  noticed  deviations  from  this  equation  in  their  critical  field 

measurements  up  to  t  ■  6.  If  the  first  term  in  eq.  (2)  is  identified  with  C  ,  it 

follows  that  eq.  (3)  is  inadequate  to  represent  the  data,  as  can  be  seen  in  fig.  1. 

In  this  graph  is  plotted  In  (C  /y  T  )  versus  T  /T.  In  the  same  graph  the  function 
3  ®  s  c  c 

0.144  T  /y  T£  is  also  drawn,  which  represents  the  data  within  the  accuracy  for  this 
term  (10  -  20$). 

4 

Cooper  proposed  that  anisotropy  in  the  energy  gap  would  lead  to  two  exponential 
temperature  dependences  for  C  .  However,  the  result  of  our  measurements  is  very 

CS 

different  from  that  expressed  in  eq.  (3).  It  does  not  appear  possible  to  represent 
the  data  with  two  exponentials,  within  reasonable  limits  for  anisotropy. 
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Ginsberg  at  al observed.  In  lead,  precursors  In  the  absorption  of  electromag¬ 
netic  waves  with  energies  about  4th  of  the  energy  gap.  One  could  assume  a  narrow  set 

3  3 

of  energy  levels  In  the  energy  gap,  but  for  such  a  set  to  lead  to  a  T  contribution 
In  the  specific  heat  appears  rather  unlikely. 

3 

A  curious  result  can  be  obtained  in  the  following  way.  A  T  -dependence  for  C 

es 

leads  thermodynamically  to  the  parabolic  law  for  the  critical  field: 

H  -  H  (1  -  (T/T  )2)  .  (4) 

o  c 

Using  the  relation  between  AC  (AC  ■  C  -  C  *  0.144  «  2,98  T)  and  d^H^/dT*  we  find 

s  n 

Ho  ■  806  Gauss,  In  excellent  agreement  with  the  value  found  from  magnetic  data 

(H  ■  802,6  Gauss).  However,  this  analysis  also  yields  a  value  of  T  ■  7.88°K,  whereas 

O  q  C 

the  normal  transition  temperature  for  lead  Is  7.19  K.  Eq.  (4)  agrees  reasonably  well 
with  magnetic  measurements  below  2°K,  but  above  this  temperature  it  leads  to  ever 
Increasing  discrepancies.  Therefore,  the  specific  heat  should  contain  other  terms 
than  the  T  one  which  appear,  however,  to  be  small.  If  one  lets  in  eq.  (4)  be 
temperature  dependent. 


Tc  -  7.88  -  In  (1  +  T/7.19)  ,  (5) 

the  magnetic  data  below  7°K  are  represented  within  1  Gauss.  This  formula  cannot  be 

2 

correct,  for  It  predicts  a  negative  T  term  In  the  specific  heat,  which  Is  not  observed 
below  1°K,  Still  the  success  of  this  temperature  dependence  of  T  ,  to  which  the 
energy  gap  Is  related,  suggests  two  possibilities: 

1.  The  lowering  of  the  energy  by  pairing  the  electrons  is  temperature  dependent  or 

2.  The  temperature  dependence  of  the  energy  gap  Is  stronger  for  lead  than  that 
calculated  by  the  BCS  theory. 

Later  measurements  of  a  single  crystal  sample  of  lead  gave  values  for  the  specific 
heat  In  the  normal  and  superconducting  state  of: 

C  -  3.03  T  +  1.61  T3  +  0.023  T5 
n 

Cs  -  0.144  T3  +  1.61  T3  +  0.028  T5 

1.  J.  Bardeen,  L.  N.  Cooper  and  J.  R.  Schrleffer,  Phys.  Rev.  I0£  (1957)  1175. 

2.  D.  L.  Waldorf  and  G.  A.  Alers,  J.  Appl.  Phys.  33  (1962)  3266. 
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3.  D.  L.  Decker,  D.  E.  Mapother  and  R.  W.  Shaw,  Phys.  Rev.  M2  (1958)  1888. 

4.  L.  N.  Cooper,  Phys.  Rev.  Letters  2  (1959)  17. 

5.  D.  N.  Ginsberg,  P.  L.  Richards  and  H.  Tlnkham,  Phys.  Rev.  Letters  _3  (1959)  337. 


THERMAL  CONDUCTIVITY  AND  THERMOELECTRIC  POWER 
N.  Pearl  man 

The  large  increase  in  thermal  resistance  which  has  been  observed  at  low  tenper- 
atures  in  Ge  doped  with  Sb  and  As  and  its  diminution  on  compensation  with  Ga  has  been 
shown  in  previous  Progress  Reports  to  depend  on  the  concentration  of  neutral  donors, 

N^  in  the  material.  The  effect  is  larger  for  Sb  donors  than  for  As  in  concentration 
ranges  where  the  Inpurity  levels  are  in  the  energy  gap  and  not  in  the  conduction  band. 
Recently,  Carruthers  and  Griffin'  have  proposed  a  mechanism  for  this  effect  which 
differs  in  some  respects  from  that  proposed  earlier  by  Keyes,  although  the  results  of 
their  calculations  are  in  general  agreement  with  those  of  Keyes  and  with  the  measure¬ 
ments.  Their  calculations  contain  predictions  which  are  of  interest  with  regard  to 
the  projected  extension  of  our  measurements  below  1°K  and  also  with  respect  to  the 
conclusions  that  can  be  drawn  as  to  the  nature  of  the  electron  levels  In  compensated 
material  with  large  impurity  concentrations. 

Keyes  calculated  the  scattering  of  phonons  as  a  result  of  the  virtual  transitions 
between  the  singlet  ground  state  and  the  triplet  excited  state,  higher  by  4A,  in  the 
neutral  donor  atoms,  using  the  elastic  strain  due  to  the  phonons  as  the  perturbation 
responsible  for  the  transitions.  The  coupling  constant  for  static  strains  is  related 
to  4A,  and  the  difference  in  this  quantity  between  Sb  and  As  accounted  for  the 
observed  difference  in  the  magnitude  of  added  thermal  resistance  for  the  two  impuri¬ 
ties.  Carruthers  and  Griffin  object  to  this  procedure  on  the  ground  that  It  is 
applicable  to  the  dynamic  phonon  strain  only  for  phonon  frequency  m  much  less  than 
46/h,  whereas  the  effect  occurs  at  temperatures  for  which  the  dominant  phonons  have 
Pm  4A/1Y.  They  therefore  calculate  the  cross-section  for  resonance  fluorescent  scat¬ 
tering  of  phonons,  involving  the  triplet  and  singlet  levels.  When  they  take  boundary 
scattering  into  account  as  well  they  find  excellent  agreement  with  the  measured  values 

between  4°K  and  l°K,  including  the  difference  observed  in  temperature  dependence  of 

4  3 

thermal  conductivity  (T  for  Sb;  T  for  As). 
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At  higher  temperatures  the  calculations  give  values  of  thermal  conductivity 
which  are  too  high,  presumably  because  of  the  neglect  of  Umklapp  scattering.  Below 
1°K  they  predict  a  rather  sharp  decrease  In  the  extra  resistance9  which  should  be 
easily  observed  within  the  planned  precision  of  our  projected  measurements. 

Although  Carruthers  and  Griffin  confine  themselves  to  the  non-compensated  case, 
their  calculation  does  contain  information  about  compensated  material  as  well.  It  Is 
based  on  the  wave-function  of  the  bound  donor  electron,  and  In  particular  the  ampli¬ 
tude  at  small  distance  from  the  core  determines  the  cut-off  in  the  extra  thermal 
resistance  at  higher  phonon  frequencies,  that  is,  at  higher  temperatures.  Now  the 
effect  In  compensated  material  has  been  shown  to  be  the  same  as  that  in  non-compen- 
sated  material  with  the  same  concentration  N^,  even  though  the  total  impurity  concen¬ 
tration  Nj  is  100  times  larger  in  the  compensated  case.  This  means  that  the  wave- 
function  does  not  differ  enough  in  the  two  cases  to  make  a  difference  in  the  calcula¬ 
tion  of  thermal  resistance.  We  have  previously  presented  evidence  from  electrical 
resistance  and  Hall  effect  measurements  that  whereas  the  donor  levels  overlap  the 
conduction  band  for  sufficiently  high  In  single-doped  Ge  they  must  be  in  the  energy 
gap  in  compensated  material  with  comparable  but  much  lower  N^.  The  calculations 
of  Carruthers  and  Griffin  present  the  possibility  that  this  earlier  qualitative  conclu¬ 
sion  may  lead  to  a  quantitative  upper  limit  on  the  permissible  difference  between  the 
compensated  Ge  with  conparable  concentrations  N^. 
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IV.  GENERAL 


X-RAY  LATTICE  PARAMETER 
R.  C.  Buschert 

Preliminary  measurements ,  to  test  the  equipment  and  the  method,  on  the  change  In 
lattice  parameter  of  germanium  due  to  doping  with  arsenic  have  been  made.  In  the 
method  used,  the  change  In  lattice  parameter  of  the  single  crystal  being  measured  is 
determined  by  measuring  the  angular  shift  in  a  diffraction  line  position  (12,  0,  0) 

between  the  sample  and  a  standard  undoped  single  dislocation  free  germanium  crystal. 

19 

The  first  sample  was  arsenic  doped  to  10  /cc  and  showed  an  angular  line  shift  of 
20"  of  arc.  This  indicates  a  lattice  parameter  increase  of  A  d/d  ■  +  5  x  10  . 

The  half  width  of  the  diffraction  line  increases  from  0.6"  of  arc  for  the  undoped 
dislocation  free  reference  crystal  to  30"  of  arc  for  the  arsenic  doped  crystal  (with 
many  dislocations).  It  is  presumed  that  most  of  this  broadening  is  due  to  disloca¬ 
tions.  However  the  possibility  exists  that  the  lattice  parameter  change  may  be 
partly  due  to  the  dislocations.  At  the  present  an  undoped  sample  with  dislocations  is 
being  prepared  In  order  to  determine  the  influence  of  dislocations  on  the  lattice 
parameter . 

In  order  to  accurately  determine  the  concentration  of  the  impurities  gallium  and 
arsenic  In  germanium  work  has  been  done  on  neutron  activation  of  these  elements  in 
the  Purdue  nuclear  reactor. 

This  bombardment  and  activity  measurement  Indicated  that  it  should  be  possible 
to  determine  the  concentration  of  either  or  both  (as  In  a  conpensated  sample)  of 
these  inpurities  In  germanium  to  a  precision  of  (1  -  5)4  for  concentrations  of 
lo'Vcc  and  greater. 


GROWTH  OF  SINGLE  CRYSTALS  OF  ZINC  TELLUR  IDE 
P.  Klose  and  H.  Y.  Fan 

Single  Crystals  of  Zinc  Tellurlde  have  been  grown  from  the  vapor  phase  and  by  the 
Bridgman  technique.  Both  methods  have  been  used  with  equal  results.  Since  the 
Bridgman  technique  permits  htgher  growth  rates  than  the  vapor  phase  amthod,  most 
crystals  have  been  grown  in  this  fashion. 
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The  crystals  have  been  checked  with  x-rays  and  optical  methods  to  determine  the 

orientation  and  size  of  the  individual  grains  of  an  ingot.  Usually  from  two  to  four 

3 

grains  made  up  an  ingot  of  about  three  cm  .  Experimental  specimens  have  been  cut  so 
as  to  avoid  grain  boundaries. 

A  mullite  pressure  container  was  used  to  contain  the  vacuum  sealed  quartz 

ampoule  which  served  as  crucible.  This  container  became  especially  useful  during  the 

doping  experiments  with  arsenic,  indium,  and  gallium  where  vapor  pressures  of  several 

atmospheres  had  to  be  sustained  by  the  quartz  anpoule. 

Measurements  of  Hall  coefficient  and  resistivity  have  been  made  on  several 

sanples  In  the  temperature  range  from  300°K  to  78°K.  The  results  are  tabulated  below. 

Figures  (26)  and  (27)  show  the  temperature  dependence  of  the  Hall  coefficient  and  the 

resistivity  of  Zinc  Telluride  sanples  doped  with  various  concentrations  of  arsenic 

163  183 

ranging  from  2.6  x  10  /cm  to  1.6  x  10  /cm  .  The  activation  energy  as  calculated 

from  the  Hall  data  varies  in  that  range  from  zero  to  0.098  eV.  Most  of  the  samples 

grown  without  doping  show  activation  energies  around  0.15  eV,  The  activation  energy 

of  copper  inpurity  has  been  reported  by  M*  Aven  to  be  0.15  eV. 

Doping  with  gallium  did  not  change  the  type  of  the  usually  p-type  crystals  but 

2  3  5  3 

changed  the  Hall  coefficient  from  10  cm  /coul  to  10  cm  /coul .  Doping  with  indium 

6  3 

produced  n-type  material  with  a  Hall  coefficient  of  10  cm  /coul.  So  far  we  have  not 
successfully  grown  n^type  crystals  of  higher  carrier  concentration.  Attempts  are 
sti 1 1  being  made. 
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Samples 

300 

°K 

78°K 

Remarks 

3 

3 

d  cm 
coul 

p  (O-cm) 

R 

coul 

p  (fl-cm) 

ZnTe  -  37 

+  100 

1.2 

+  6.9  X  105 

1 .6  x  105 

Vapor  phase-As 

ZnTe  -  38 

+  425 

4.94 

Vapor  phase 

ZnTe  -  59 

+  77.5 

0.67 

+  1 .9  x  108 

7.4  x  104 

Bridgman 

ZnTe  -  111 

+  285 

2.5 

Bridgman 

ZnTe  -  125 

+  217 

2.5 

+  1.2  x  109 

Bridgman 

ZnTe  -  133 

+  84.4 

1.0 

+  3.6  x  107 

9.2  x  104 

Bridgman 

ZnTe  -  141 

+  1140 

24.1 

Bridgman 

ZnTe  -  145 

+  40.8 

0.51 

+  2.6  x  104 

65.3 

Brldgman-As 

ZnTe  -  146 

+  196 

1.54 

Vapor  phase 

ZnTe  -  147 

+  2.93  x  105 

6.2  x  I04 

Brldgman^Ga 

ZnTe  -  150 

-  1 .3  x  106 

2.48  x  104 

Br Idgman-ln 

ZnTe  -  151 

-  2.9  x  106 

3.86  x  104 

Br fdgman-ln 

ZnTe  -  152 

+  244 

2.68 

+  5.3  x  108 

Brldgman-As 

ZnTe  -  154 

+  4.1 

0.16 

+  3.3 

0.47 

Brldgman-As 

ZnTe  -  155 

+  55.5 

0.66 

+  2.64  x  105 

8.95  x  103 

Brldgman-As 
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CAPTIONS  FOR  FIGURES  20-22 


Fig.  20.  EPR  spectrum  of  Center  (II,  III)  with  H  ||  <  110  >.  Sample  A  (Vacuum 

Floating  zone  silicon,  0  ■  1 .8  x  10*®  nvt  at  ORNL)  annealed  at  100°C  for 
30  minutes.  T  m  300°K,  v  ■  24.552  108  kMc.  The  principal  lines  are  marked 
according  to  whether  they  belong  to  Set  II  or  Set  III.  The  hyperflne  satel¬ 
lites  (5.1ft  relative  amplitudes)  due  to  the  Si2?  nuclei  In  the  nearest 
lattice  sites  and  the  "triplet”  structure  (23.6ff  relative  amplitudes)  due  to 
those  In  next  nearest  sites  are  visible  on  one  of  the  Isolated  branches. 


Fig.  21.  Hyperflne  splitting  of  Center  (II,  III).  The  separation  between  the  hfs 

satellites  (5. Iff  relative  amplitude)  and  the  corresponding  branches  of  the 
fine  structure  Is  plotted  vs  the  orientation  of  the  external  field  In  a 
(110)  plane.  Points  taken  from  the  ms  ■  0  to  +  I  transitions  fall  on  the 
same  curves  and  are  not  distinguished. 


Fig.  22.  Angular  variation  plot  of  Center  (II,  III)  spectrum.  Sanple  2  (Crucible 
grown  silicon,  0  m  1.4  x  10*®  nvt  at  ORNL),  unannealed.  T  ■  300°K, 
v  ■  24.112  101  kMc.  The  external  magnetic  field,  shown  In  proton  resonance 
frequency  (Me),  Is  rotated  in  a  (110)  plane.  The  open  circles  belong  to 
Set  II  and  the  closed  circles  to  Set  III.  The  diameter  of  the  circles 
corresponds  to  the  full  width  between  inflections  of  the  resonance  lines. 

The  thick  lines  represent  the  fine  structure  lines  and  the  thin  lines  the 
hyperflne  satellite  lines.  The  branches  are  coded  by  solid,  dashed,  and 
dotted  lines  to  indicate  the  pairs  belonging  to  the  same  orientation  of  the 
center.  The  "triplet”  lines  are  also  shown  on  two  upper  branches  of  Set  II. 


FIGURE 


FIGURE 


war 


FIGURE  3 


FIGURE 


FIGURE 


T 


T 


CL 

a 


800  v/cm 


FIGURE  6 


FIGURE  7 


EQUIVALENT  TEMPERATURE 


10 


E  (volts/cm) 
FIGURE  8 


MOBILITY  (cm*/volt  sec) 


E  (volt/cm) 


FIGURE  9 


FIGURE  10 


-GaSb(Te) 
T  =  4.2°K 


o 


r 


R  (cm’cour1) 


7 T  (lOT'cmVd) 


FIGURE  17 


|Rh|  (cm* /coul) 


FIGURE  18 


FIGURE  19 


FIGURE  20 


AH  (PROTON  RES.  FREQ.,  mc/s) 


GAUSS 


V  aidues  *  (kuio)  u*dw»03  uoutUoiqv 


1 


FIGURE  23 


50, 


Vibrational  Spectrum  of 
Oxygen  in  Germanium 


T~I0-I5°K 


850 


_ I _ 

860 

Wavenumber 


(cm f) 


870 


FIGURE  24 


OlK 


Resistivity  (  Ohm  -  cm ) 


FIGURE  27 


3urdua  Hi l varsity 
DA  31-12 >~AS0(D)~|7 


Semiannual  Report 
1  October,  1962  -  31  March,  1963 


D I  ree  :or 

National  Bureau  of  Standard* 
t3oul d^r,  Colorado 
Attn:  Radio  Library 

Commander,  Air  Force  Command  and  Control  Development  Division 

Air  ^search  and  Development  Command 

iJn  So  Air  Force 

L.  G„  Hanscom  Field 

Bedford,  Massachusetts 

Attn:  CA07R  26  I 

;)« reccor 

'J  Sn  Na</al  Research  Laboratory 
Washington  25,  D*  C0 

Attn:  Csde  2027  1 

Chief  Signal  Officer 
department  of  the  Army 
Washington  25,  0o  C. 

Attn:  SIGRD  l 


No,  of  Copies 

l 


Commander 

Armed  Services  Technical  Information  Agency 
Arlington  Hall  Station 

Arlington  12,  Virginia  10 

Commander 

Wright  Air  Development  Division 
/irlght-Patterson  Air  Force  Base,  Ohio 
Attn:  WC0S1-3 
Attn:  WCLKT-4 

Department  of  the  Navy 

Bureau  of  Ships 
Sfjmlconcluctor  Unit,  Code  691A 
Washington  25,  D0  C„ 

Attn:  Nr,  A0  Hc  Young 

Cnisf,  Bureau  of  Ships 
Department  of  the  Navy 
Washington  25,  D„  C. 

Attn:  Code  327 

Office  of  Naval  Research 
Department  of  the  Navy 
Washington  25,  D„  C« 
tltn:  Code  427 


I 

I 


* 

% 

DA  3M24-AR0(D)-I7 

rt9.  fif  faftlM 

Director 

U.  S.  Naval  Research  Laboratory 

Washington  25,  D.  C. 

Attn:  Code  3560 

1 

Advisory  Group  on  Electron  Devices 

346  Broadway 

New  York  13,  New  York 

2 

Commanding  Officer 

Watertown  Arsenal 

Watertown,  Massachusetts 

Attn:  OHIO 

1 

Commander 

Air  Research  and  Development  Command 

Andrews  Air  Force  Base 

Washington  25,  D.  C. 

Attn:  RDRR-2 

1 

Commanding  Officer 

U.  S.  Army  Signal  Electronics  Research  Unit  * 

Mountain  View,  California 

1 

Commandant 

Anti-aircraft  Artillery  and  Guided  Missile  School 

Fort  Bliss,  Texas 

1 

U.  S.  Naval  Inspector  of  Ordnance 

Applied  Physics  Laboratory 

Johns  Hopkins  University 

Silver  Spring,  Maryland 

1 

Donald  K.  Stevens 

Asst.  Director  for  Metallurgy  and  Materials  Programs 

Division  of  Research 

U.  S.  Atomic  Energy  Commission 

Washington  25,  D.  C. 

1 

Hew  York  Naval  Shipyard 

Brooklyn  1,  New  York 

Attn:  Material  Laboratory,  Code  912b 

1 

Naval  Ordnance  Laboratory 

Silver  Spring,  Maryland 

Attn:  Solid  State  Division 

1 

Director,  U.  S.  Naval  Research  Laboratory 

Washington  25,  D.  C. 

Attn:  Code  6451 

1 

DA  31-114-AR0(D)-17 


Coinroanding  Officer  and  Director 
U,  S#  Naval  Electronics  Laboratory 
San  Die$o  52,  California 

Technical  Library 
OASD  (R  and  E) 

Rm  3E1065,  The  Pentagon 
Washington  25,  D.  C„ 

Chief  of  Ordnance 
Washington  25,  D.  C* 

Attn:  ORDTX-AR 

Contending  General 
U,  S.  Army  Ordnance  Missile  Command 
Redstone  Arsenal,  Alabama 
Attn:  0RDXFM5R 

Commanding  Officer 
Frankford  Arsenal 
Philadelphia  37,  Pennsylvania 
Attn:  ORDBA-FEL 

RADC  ( RAALD ) 

Griff  I  ss  Air  Force  Base,  New  York 
Attn:  Documents  Library 

Coranandlng  Officer 
Diamond  Ordnance  Fuze  Laboratory 
Connecticut  Avenue  and  Van  Ness  Street 
Washington  25,  D,  C. 

Attn:  T0  T«  Lflmatalnen 

Commanding  Officer 

U.  S0  Army  Signal  Equipment  Support  Agency 
Fort  Monmouth,  New  Jersey 
Attn:  SIGFM/ES-ASA 

Army  Liaison  Office 

U.  S*  Naval  Research  Laboratory 

Washington  25,  D«  C. 

Commander 

Rome  Air  Development  Center 
Grlfflss  Air  Force  Base,  New  York 
Attn:  RCRWE-3 

Chief  of  Ordnance 
Washington  25,  D.  C. 

Attn:  ORDTB-Materlals 


DA  31-I24-AK0(D)-17 

Continental  Army  Command  Liaison  Office 

U.  S.  Army  Signal  Research  and  Development  Laboratory 

Fort  Monmouth,  New  Jersey 

U*  S.  Office  of  Naval  Research 
National  Bureau  of  Standards 
Corona,  California 

Conmandlng  Officer 

Ur  S«  Army  Signal  Research  and  Development  Laboratory 

Fort  Monmouth,  New  Jersey 

Attn:  Technical  information  Division 

FOR  RETRANSMITTAL  TO:  Canadian  Liaison  Officer 

Office  of  the  Chief  Signal  Officer 
Room  2B274,  Pentagon 
Washington  25,  D,  C. 

Commend ing  Officer 

Engineering  Research  and  Development  Laboratories 

Fort  Belvolr,  Virginia 

Attn:  Technical  Intel  1 Igence  Branch 

Dattel Ic  Memorial  Institute 
505  King  Avenue 
Co!urri>us  I,  Ohio 
Attn:  A*  C.  Beer 

Chicago  Midway  Laboratories 
University  of  Chicago 
Chicago  37,  Illinois 
Attn:  Librarian 

Technical  Reports  Collection 
303A  Pierce  Hall 
Harvard  University 
Cambrldra  38,  Massachusetts 
Attn:  Librarian 

Franklin  Institute 

20th  Street  and  Benjamin  Parkway 

Philadelphia,  Pennsylvania 

University  of  Illinois 
Urbana,  Illinois 
Attn:  Dr.  J.  Bardeen 

Westinghouso  Electric  Corporation 
Youngwood,  Pennsylvania 
Attn:  R.  K.  Riel 


TA  3M24-AR0(D)-17 

i:orthv;cstern  University 
Evanston.  Hit  no  Is 
..ten:  Dr*  V/*  Nuxford 

Chief  of  Research  and  Development 
ACS,  leportibant  of  the  Array 
Washington  25f  D.  C. 

Mug  he*.  Aircraft  Company 
•Semiconductor  Division 
\  0.  Box  278 
'lawport  333ch,  California 
Attn:  Library 

National  Bureau  of  Standards 
!J*  S.  Department  of  Commerce 
Washington  25,  D.  C# 

-ttn:  P.  J*  Solgin 

Phllco  Corporation 
:7C”  and  Tioga  Streets 
Philadelphia,  Pennsylvania 
Attn:  Mr.  W#  H*  Forster 

.jenaral  Electric  Company 
Semiconductor  Products  Department 
;] act  ionics  Park 
^/ractjse.  New  York 

Attn:  Dr.  H.  M.  Sullivan,  Advanced  Semiconductor  Lab. 

A&ythoon  Company 
150  California  Street 
.tev/ton,  Massachusetts 
•  ttn:  Mr.  K.  L.  Gaudreau,  Administrator 
Government  Contracts 

A?)]  Telephone  Laboratories 
Murray  Will,  iNew  Jersey 
Attn:  Mr.  A.  Anderson 

rtadlo  Corporation  of  America 
Harrison,  Now  Jersey 
Attn:  C„  M.  Bost 

irown  University 
Providence  12,  Rhode  Island 

Massachusetts  Institute  of  Technology 
■>.  0.  Box  73 

-exlngton  73 ,  Massachusetts 
5 ttn:  Library  A-229 


Wo.  of  Coola» 


I 


I 


I 


1 


I 


I 


1 

I 

I 

I 


1 


r A  3!-!?.4W».flO<D)-17 


Wo.  of  Cop»a» 


!!orthv:as:arn  University 
Ivans  I  on,  Illinois 
♦  .fin:  Dr.  W.  Nuxford 

Chief  of  ^search  and  Development 
ACS,  Eopnrtn&nt  of  the  Army 

Joshing  ton  25,  D„  C# 

(iugSies  Aircraft  Company 
oeiniconductor  Division 
0o  Bo:c  278 

Newport  3each,  California 
Attn:  Library 

National  Bureau  of  Standards 
;J.  S„  Department  of  Conmarce 
Washington  25 ,  D,  C. 

»ttn:  Pa  J.  Sotgin 

Phltco  Corporation 
"C"  and  Tioga  Streets 
’hUacblphia,  Pennsylvania 
Attn:  Mr.  V/*  H.  Forster 

oonanal  Electric  Company 
Cemlconduetor  Products  Department 
•Electronics  Park 
Syracuse,  New  York 

Attn:  Dr.  H.  M.  Sullivan,  Advanced  Semiconductor  Lab. 

Asythaon  Company 
150  California  Street 
Lawton,  Massachusetts 
\ttn:  Mr.  K.  L.  Gaudreau,  Administrator 
Government  Contracts 

Wall  Telephone  Laboratories 
Murray  Mill,  New  Jersey 
Attn:  Mr.  A.  Anderson 

rtadio  Corporation  of  America 
Harrison,  New  Jersey 
Attn:  C*  M.  Bost 


1 


1 


I 


1 


I 


I 


l 


I 


1 


drown  University 

Providence  12,  Rhode  Island  1 

Massactiusetts  Institute  of  Technology 
■\  0.  Box  73 

.exlngtcn  73,  Massachusetts 
5ttn:  Library  A-229 


1 


DA  3M24-AR0(D)  17 


No.  of  Conies 


University  of  Hlnneso 
Institute  of  Technology 
Minneapolis  14,  Minnesota 
Attn:  Dr.  A.  ven  der  Ziel 

New  York  University 
College  of  Engineering 
University  Heights 
New  York  53,  New  York 
Attn:  Dr.  I.  Cadoff 

Massachusetts  Institute  of  Technology 
Research  and  Development  Laboratories 
Cambridge  39,  Massachusetts 
Attn:  Documents  Office,  Room  2 OB-22 I 

Notre  Dame  University 
Physics  Department 
South  Bend,  Indiana 
Attn:  Prof.  J.  C.  Buch 

University  of  Michigan 
Willow  Run  Research  Center 
Ypsilanti,  Michigan 
Attn:  Librarian,  SR IA 

Pennsylvania  State  University 
College  of  Engineering  and  Architecture 
University  Park,  Pennsylvania 
Attn:  S.  H.  Chamberlain,  Librarian 

Sanders  Associates,  Inc. 

Nashua,  New  Hampshire 

North  American  Aviation 
International  Airport 
Los  Angeles,  California 

National  Carbon  Company 

Division  of  Union  Carbide  Corporation 

P.  0.  Box  61  If 

Cleveland,  Ohio 
Attn:  Library 

Mr.  Robert  K.  Willardson,  Assistant  Chief 

Physical  Chemistry  Division 

Battel le  Memorial  Institute 

505  King  Avenue 

Colucfcus  I,  Ohio 


1 


1 


1 


1 


I 

1 


I 


1 


I 


DA  3l-124-AR0(D)«t7 


Wo.  of  Copies 


Honeywell  Research  Center 
500  Washington  Avenue,  South 
Hopkins,  Minnesota 

Attn:  Dr.  Van  W.  Bearlnger  I 

Semiconductor  Components  Library 
Texas  Instruments,  Inc. 

P.  0.  Box  5012 

Dallas  22,  Texas  I 

Stanford  University 
Palo  Alto,  California 

Attn:  Dr.  John  Moll,  Electrical  Engineering  Department  I 

Commander 

U,  S.  Naval  Ordnance  Test  Station 
China  L  .ice,  California 

Attn:  Code  5019  | 

International  Rectifier  Corporation 
1521  East  Grand  Avenue 
El  Segundo,  California 

Attn:  Dr.  C.  A,  Escoffery  1 

Chief,  U.  S.  Army  Security  Agency 
Arlington  Hall  Station 

Arlington  12,  Virginia  2 

Commanding  Officer 

U.  S.  Army  Research  Office  (Durham) 

Attn:  CRD-AA-1P,  Mr.  Ulsh 
Box  CH,  Duke  Station 

Durham,  N;  Carolina  3 

Commanding  Officer 

U.  S.  Army  Signal  Research  and  Development  Laboratory 
Fort  Monmouth,  New  Jersey 

Attn:  Director  of  Research  I 

Attn:  Chief,  Technical  Documants  Center  j 

Attn:  Chief,  Technical  Information  Division  3 

Attn:  Exploratory  Research  Division  "E"  1 

Attn:  Microwave  Quantum  and  Electronics  Br,  Solid  State  and  Frequency  Control  Dlv  1 

Attn:  Mr.  P.  Newman,  Solid  State  and  Frequency  Control  Dlv  3 

Deputy  President 

U.  S,  Army  Security  Agency  Board 
Arlington  Hall  Station 
Arlington  12,  Virginia 


I 


DA  31-I24-AR0(D)-I7 


Sylvanla  Electric  Products,  Inc. 
Semiconductor  Division 
100  Sylvan  Road 
Woburn,  Massachusetts 
Attn:  Mr.  F.  H.  Bower 

Diamond  Ordnance  Fuze  Laboratories 
U,  S.  Army  Ordnance  Corps. 

Attn:  ORD  and  L  -  450-636 
Mr.  Raymond  H.  Comyn 
Washington  25,  D.  C. 

Total  number  of  copies  to  be  distributed 


This  contract  is  supervised  by  the  Institute  for  Exploratory  Research 
Scientific  contact  at  Purdue:  Dr.  H.  Y.  Fan 

Scientific  contact  at  Signal  R  and  D  Laboratory:  Dr,  S.  Benedict  Levin 

Telephone  201-535-1308 


